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Abstract:  Recent  observational  and  theoretical  studies  of  the  global  properties  of  small-scale  atmospheric  gravity  waves  have 
highlighted  the  global  effects  of  these  waves  on  the  circulation  from  the  surface  to  the  middle  atmosphere.  The  effects  of  gravity 
waves  on  the  large-scale  circulation  have  long  been  treated  via  parametrizations  in  both  climate  and  weather  forecasting  applications. 
In  these  parametrizations,  key  parameters  describe  the  global  distributions  of  gravity  wave  momentum  flux,  wavelengths,  and 
frequencies  of  the  waves.  Until  recently,  global  observations  could  not  define  the  needed  parameters  because  the  waves  are  small 
in  scale  and  intermittent  in  occurrence.  Recent  satellite  and  other  global  data  sets  with  improved  resolution  along  with  innovative 
analysis  methods  are  now  providing  constraints  for  the  parametrizations  that  can  improve  the  treatment  of  these  waves  in  climate 
prediction  models.  Research  using  very  high  resolution  global  models  has  also  recently  demonstrated  the  capability  of  resolving 
gravity  waves  and  their  circulation  effects,  and  when  tested  against  observations,  these  models  are  showing  some  very  realistic 
properties.  Here  we  review  recent  studies  on  gravity  wave  effects  in  stratosphere-resolving  climate  models,  recent  observations  and 
analysis  methods  that  reveal  global  patterns  in  gravity  wave  momentum  fluxes,  and  results  of  the  very  high  resolution  model  studies, 
and  we  outline  some  future  research  needs  to  improve  the  treatment  of  these  waves  in  climate  simulations.  Copyright  ©  2009  Royal 
Meteorological  Society 
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1  Introduction 


Small-scale  atmospheric  waves,  called  gravity  waves  (or 

‘Correspondence  to:  M.  J.  Alexander,  NWRA/Colorado  Research  Asso¬ 
ciates  Division,  3380  Mitchell  Lane,  Boulder,  CO,  80301,  USA  internal  waves),  have  sources  in  the  troposphere  such 
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as  flow  over  topography,  convection,  and  jet  imbalance. 
As  these  waves  propagate  upward  and  dissipate,  they 
force  the  atmospheric  circulation  at  altitudes  near  the 
tropopause  and  well  above,  in  the  stratosphere  and  meso¬ 
sphere.  Global  circulation  models  used  for  weather  and 
climate  forecasting  may  resolve  some  gravity  waves  but 
include  the  effects  of  unresolved  gravity  wave  forcing 
via  parametrizations.  The  parametrizations  compute  a 
momentum  forcing  term  by  making  assumptions  about 
the  unresolved  wave  properties  that  have  not  been  well- 
constrained  by  observations.  The  assumptions  are  formu¬ 
lated  as  a  set  of  tuning  parameters  that  are  used  to  adjust 
the  circulation  and  temperature  structure  in  the  upper  tro¬ 
posphere  and  middle  atmosphere. 

Different  models  have  different  goals  and  apply  dif¬ 
ferent  types  of  parametrizations  with  unique  tunings. 
For  example,  chemistry-climate  models  use  gravity  wave 
parametrizations  to  adjust  their  stratospheric  circula¬ 
tion  and  polar  temperatures  to  realistic  values,  impor¬ 
tant  for  accurate  simulation  of  ozone  chemistry.  Climate 
and  weather  forecasting  models  have  traditionally  used 
mountain  wave  parametrizations  to  improve  the  structure 
of  the  winter  jets  and  horizontal  temperature  gradients 
near  the  tropopause,  and  improve  surface  wind  distribu¬ 
tions.  Recent  studies  have  highlighted  additional  impor¬ 
tant  effects  of  gravity  waves  in  climate  models  that  under¬ 
score  the  need  to  reduce  uncertainties  in  gravity  wave 
parametrizations.  The  requirements  needed  from  observa¬ 
tions  to  constrain  the  tuning  parameters  in  gravity  wave 
parametrization  schemes  have  been  impossible  to  achieve 
on  the  global  scale  until  fairly  recently. 

What  is  needed  is  characterization  of  the  vertical  flux 
of  horizontal  pseudomomentum,  and  three-dimensional 


(3D)  wave  propagation  properties.  We  will  use  the  com¬ 
mon  approximation  for  pseudomomentum  flux  as  momen¬ 
tum  flux,  which  for  gravity  waves  can  be  written  as  the 
product  of  density  and  the  covariance  of  horizontal  u' 
and  vertical  w'  wind  perturbations.  It  is  often  written  as 
p{u'w' ,  v'w'),  the  overbar  represents  a  spatial  or  temporal 
average  over  a  wavelength  or  period.  The  3D  wave  proper¬ 
ties  include  horizontal  and  vertical  wavelength,  as  well  as 
propagation  direction.  Only  in  the  last  decade  have  high- 
resolution  satellite  observations  allowed  the  estimation  of 
these  wave  properties  and  momentum  fluxes.  However, 
each  individual  satellite  measurement  technique  can  only 
observe  a  portion  of  the  full  3D  spectrum  of  gravity  waves. 
Errors  in  the  momentum  fluxes  derived  from  global  obser¬ 
vations  remain  large,  although  local  case  studies  can  now 
be  quite  accurate.  Some  very  high  resolution  global  mod¬ 
elling  studies  that  include  a  middle  atmosphere  are  now 
resolving  much  more  of  the  gravity  wave  spectrum  and 
are  being  used  to  study  the  wave  sources,  propagation 
and  dissipation,  and  the  momentum  forcing  of  the  circula¬ 
tion  without  the  use  of  any  gravity  wave  parametrizations. 
In  recent  years,  new  methods  for  inferring  the  missing 
momentum  forcing  due  to  unresolved  waves  in  lower  res¬ 
olution  models  have  also  been  developed  using  advanced 
data  assimilation  methods. 


With  these  recent  developments,  we  can  collectively 
better  constrain  and  quantify  the  role  of  gravity  waves 
in  the  momentum  budget  of  the  global  circulation.  The 
purpose  of  this  review  is  to  provide  an  overview  of  recent 
developments  in  this  rapidly  developing  field. 
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2  Gravity  wave  issues  in  stratosphere-resolving 
global  climate  models 

General  circulation  models  used  for  climate  studies  of  the 
troposphere-stratosphere  system  have  generally  used  rel¬ 
atively  coarse  resolution  (~  2-5°  in  the  horizontal  and 
variable  in  the  vertical  but  usually  ~  3km  in  the  strato¬ 
sphere)  and  prescribed  sea  surface  temperatures  (SSTs). 
Depending  upon  the  problem  at  hand,  they  may  be  either 
dynamics-only  models  (e.g.,  those  that  use  prescribed 
ozone  fields  to  compute  solar  heating)  or  models  coupled 
to  a  fully  interactive  chemistry  scheme  (e.g.,  those  that 
use  predicted  ozone  to  compute  solar  heating).  Given  that 
multi-year  simulations  are  necessary,  their  coarse  resolu¬ 
tion  is  dictated  by  current  computional  resources,  a  situ¬ 
ation  that  is  not  likely  to  change  much  in  the  foreseeable 
future. 

While  their  coarse  resolution  is  generally  fine  enough 
to  capture  most  of  the  Rossby  wave  activity  in  the  atmos¬ 
phere,  it  is  not  fine  enough  to  resolve  the  bulk  of  the 
small-scale  gravity  waves  that  play  a  large  role  in  the 
momentum  budget  of  the  middle  atmosphere.  To  account 
for  the  effects  of  these  unresolved  disturbances,  gravity 
wave  drag  (GWD)  parameterzations  must  be  used.  Fol¬ 
lowing  the  pioneering  work  of  Lindzen  (1981)  and  sub¬ 
sequent  refinements  by  Holton  (1982),  the  first  imple¬ 
mentation  of  GWD  parametrizations  in  global  climate 
models  was  in  the  context  of  orographically-forced  grav¬ 
ity  waves  (Palmer  et  ah,  1986;  McFarlane,  1987).  How¬ 
ever,  as  model  lid  heights  increased  into  the  mesosphere, 
parametrizations  for  non-orographic  gravity  waves  (e.g., 
gravity  waves  forced  by  convection  and  other  sources) 
became  necessary.  As  a  result,  GWD  parametrizations 
are  generally  grouped  into  two  categories:  1)  orographic 
GWD  (OGWD)  schemes,  which  represent  zero  phase 
speed  waves  forced  by  subgridscale  topography,  and  2) 


non-orographic  GWD  (NGWD)  schemes,  which  represent 
non-zero  phase  speed  waves  forced  by  mechanisms  other 
than  flow  over  topography.  See  McLandress  (1998)  for  a 
review  of  a  number  of  those  parametrizations  and  their 
impact  in  climate  model  simulations. 

In  order  to  make  the  GWD  parametrization  prob¬ 
lem  tractable  and  computationally  efficient,  a  number 
of  important  simplifying  assumptions  must  be  made 
about  the  waves  and  their  sources,  including  the  spec¬ 
trum  of  momentum  flux.  Assumptions  common  to  current 
parametrizations  are  that  the  gravity  waves  propagate  only 
vertically  and  instantly  through  the  column,  and  that  the 
waves  propagate  in  a  conservative  fashion  up  to  a  height 
at  which  an  amplitude  threshold  is  exceeded,  the  wave 
breaks,  deposits  momentum  flux,  and  generates  GWD  (see 
references  in  McLandress,  1998).  All  GWD  parametriza¬ 
tions  employ  tunable  parameters  which  act  to  scale  the 
wave  drag  and/or  change  the  breaking  heights.  In  princi¬ 
ple,  these  parameters  should  be  constrained  by  observa¬ 
tions.  In  practice,  however,  they  are  determined  by  tuning 
the  models  to  observed  climatologies. 

In  the  following  three  subsections,  we  discuss  a  num¬ 
ber  of  recent  issues  concerning  the  impact  of  parameter¬ 
ized  GWD  on  global  simulations  of  climate  and  climate 
change.  We  focus  on  the  region  below  the  stratopause 
since  that  is  where  these  recent  advances  have  occurred. 

2. 1  Issues  regarding  angular  momentum  conservation 

The  transfer  of  momentum  by  gravity  waves  is  an  impor¬ 
tant  process  in  the  global  angular  momentum  budget  of 
the  atmosphere.  Parametrizations  of  this  transfer  must 
respect  the  principle  of  conservation  of  angular  momen¬ 
tum,  which  represents  an  important  constraint  on  the 
dynamics.  Since  current  GWD  parametrizations  assume 
that  momentum  transfers  are  strictly  vertical,  angular 
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momentum  conservation  is  equivalent  to  zonal  momen¬ 
tum  conservation.  Shepherd  and  Shaw  (2004)  formulated 
a  momentum  conservation  constraint  for  GWD  parame¬ 
terization  in  a  model  vertical  column,  which  requires  that 
any  gravity  wave  momentum  flux  across  a  vertical  model 
level  must  be  entirely  absorbed  in  the  atmosphere  above. 
The  constraint  expresses  the  fact  that  the  momentum  bud¬ 
get  is  closed  with  regard  to  the  upper  boundary  (i.e.,  there 
can  be  no  radiation  of  gravity  wave  momentum  flux  to 
space). 

The  implications  of  the  gravity  wave  momentum 
conservation  constraint  for  modelled  climate  have  been 
explored  by  Shaw  and  Shepherd  (2007)  and  Shaw  et  al. 
(2009).  Using  a  zonally-symmetric  mechanistic  model, 
Shaw  and  Shepherd  (2007)  showed  that  violating  momen¬ 
tum  conservation  leads  to  a  non-robustness  in  the  response 
of  parametrized  nonorographic  gravity  waves  to  a  radia¬ 
tive  perturbation  in  the  middle  atmosphere.  In  particu¬ 
lar,  when  momentum  was  not  conserved,  the  response 
became  sensitive  to  changes  in  the  gravity  wave  source 
spectrum,  parametrization  scheme,  background  flow  and 
model  lid  height.  Momentum  conservation  was  violated 
by  either  letting  gravity  wave  momentum  flux  escape  to 
space  or  through  the  use  of  a  Rayleigh  friction  layer.  The 
sensitivity  to  model  lid  height  resulting  from  nonconser¬ 
vation  explained  the  sensitivity  to  lid  height  found  by 
Lawrence  (1997)  and  is  unacceptable  for  model  intercom¬ 
parison  purposes.  The  spurious  response  resulting  from 
nonconservation  was  also  non-negligible  in  terms  of  its 
magnitude  and  downward  influence  into  the  troposphere. 
The  downward  influence  is  described  by  the  “downward 
control”  mechanism  (Haynes  et  al.,  1991),  whereby  wave 
drag  in  the  stratosphere  induces  meridional  circulation 
cells  that  extend  into  the  troposphere  and  modify  tempera¬ 
ture  gradients,  zonal  wind  patterns,  and  sea-level  pressure. 


When  momentum  was  conserved  the  response  became 
robust  to  changes  in  the  above  mentioned  factors. 

Using  a  comprehensive  global  climate  model,  Shaw 
et  al.  (2009)  quantified  the  impact  of  conserving  both  oro¬ 
graphic  and  nonorographic  gravity  wave  momentum  flux 
on  3D  modelled  climate  and  on  the  climate  response  to 
idealized  ozone  depletion.  When  the  model  lid  was  placed 
at  10  hPa,  nonconservation  led  to  large  biases  in  the  mean 
climate  and  its  variability,  which  extended  to  the  sur¬ 
face.  The  errors  due  to  nonconservation  resulted  not  only 
from  the  direct  effect  associated  with  the  missing  force 
attributed  to  the  missing  gravity  wave  momentum  flux  but 
also  from  errors  in  the  resolved  planetary  waves  in  the 
model.  The  effects  on  the  planetary  waves  resulted  from  a 
feedback  through  the  modelled  winds,  which  were  always 
closer  to  reality  when  momentum  was  conserved.  In  prac¬ 
tice  the  errors  which  arise  from  nonconservation  could  be 
corrected  by  re-tuning  the  gravity  wave  parametrization 
parameters,  however  this  would  be  tuning  against  an  error. 
When  the  model  was  perturbed  by  idealized  ozone  deple¬ 
tion  there  were  large  biases,  both  in  the  magnitude  and 
in  the  sign  of  the  planetary  wave  and  parametrized  grav¬ 
ity  wave  responses  relative  to  previous  modelling  studies 
(Manzini  et  al.,  2003).  When  the  model  lid  was  placed  at 
0.001  hPa,  there  was  little  impact  of  momentum  conser¬ 
vation  on  both  the  mean  climate  and  its  response  to  the 
idealized  ozone  depletion,  which  is  expected  due  to  the 
small  amount  of  gravity  wave  momentum  flux  reaching 
0.001  hPa. 

Momentum  conservation  in  GWD  parametrization  is 
an  important  theoretical  constraint.  It  is  most  useful  when 
the  total  amount  of  gravity  wave  momentum  flux  input 
into  the  parametrization  is  accurate.  Observations  of  grav¬ 
ity  wave  momentum  flux  are  key  to  constraining  GWD 
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parametrizations  and  in  conjunction  with  momentum  con¬ 
servation  ensure  the  physical  realism  of  the  large-scale 
response  to  the  GWD. 

2.2  Issues  regarding  orographic  GWD 

The  need  to  parametrize  the  effects  of  small-scale  grav¬ 
ity  waves  generated  by  unresolved  topography  in  global 
models  has  been  acknowledged  for  over  two  decades  (e.g.. 
Palmer  et  ah,  1986;  McFarlane,  1987).  Inclusion  of  their 
effects  leads  to  substantial  improvements  in  simulations 
of  the  mid-latitude  lower  stratosphere  in  northern  winter, 
and,  by  remote  influence,  the  troposphere.  This  occurs  as 
follows:  The  bulk  of  the  OGWD  occurs  in  the  mid-latitude 
lower  stratosphere,  where  it  leads  to  a  deceleration  of  the 
zonal  wind  in  that  region.  These  weakened  winds  cause  a 
separation  of  the  subtropical  and  polar  night  jets,  which 
in  turn  alters  the  refractive  properties  for  planetary  waves, 
leading  to  changes  in  the  patterns  of  resolved  wave  drag 
(Sigmond  and  Scinocca,  2010).  The  impact  of  both  the 
parametrized  and  resolved  wave  drag  on  the  troposphere 
then  follows  from  “downward  control”  (Haynes  et  ah, 
1991).  A  number  of  recent  studies  indicate  that  the  “down¬ 
ward  control”  signal  in  the  troposphere  can  be  amplified 
by  feedback  mechanisms  involving  tropospheric  eddies 
(e.g.  Song  and  Robinson,  2004;  Chen  and  Zorita-Gotor, 
2008). 

Thus,  in  order  to  get  the  stratospheric  and  tropo¬ 
spheric  responses  to  OGWD  right  for  the  right  reas¬ 
ons,  it  is  important  that  the  tunable  parameters  be  con¬ 
strained  by  observations.  A  common  tunable  parameter 
in  OGWD  schemes  is  a  so-called  efficiency  factor  (e), 
which  acts  to  scale  the  drag  without  affecting  break¬ 
ing  heights.  The  Whole  Atmosphere  Community  Climate 
Model  (WACCM;  Garcia  et  al.  2007)  uses  a  variant  of 
the  McFarlane  (1987)  scheme,  and  calls  the  efficiency 
factor  e.  In  the  Scinocca  and  McFarlane  (2000)  scheme, 


the  efficiency  factor  is  included  in  the  parameter  G.  The 
values  of  these  parameters  are  typically  chosen  to  be 
those  which  produce  the  smallest  biases  in  quantities  such 
as  zonal  mean  zonal  winds  and  mean  sea-level  pressure 
(e.g..  Palmer  et  ah,  1986,  McFarlane,  1987,  Scinocca  et 
al.,  2008).  Separate  parameters  can  change  the  breaking 
heights  (e.g.,  “Frcrjt  in  Scinocca  and  McFarlane,  2000). 
Observations  are  not  likely  to  be  able  to  directly  con¬ 
strain  such  parameters,  but  can  constrain  physical  quan¬ 
tities  such  as  momentum  flux  and  intermittency.  (See  Sec¬ 
tion  3.)  From  parametrizations  in  models,  it  is  possible  to 
reconstruct  and  output  these  physical  quantities  so  that  in 
future  studies,  the  tuning  parameters  can  be  constrained 
by  observations  through  this  indirect  route. 

Figure  1  shows  climatological  zonal  mean  zonal 
winds  for  December-February  from  the  Canadian  Centre 
for  Climate  Modelling  and  Analysis  (CCCma)  model 
(which  uses  the  Scinocca  and  McFarlane  (2000)  scheme), 
computed  using  two  values  of  G  (Sigmond  and  Scinocca, 
2010).  The  “weak”  setting  (left)  produces  more  realistic 
zonal  mean  winds  in  the  lower  stratosphere,  but  at  the 
expense  of  less  realistic  mean  sea-level  pressures  (not 
shown).  The  “strong”  setting  (middle)  produces  overly 
weak  winds  in  the  lower  stratosphere,  but  more  realistic 
mean  sea-level  pressures.  The  wind  differences  (right)  are, 
in  fact,  accentuated  when  the  model’s  upper  boundary  is 
raised  to  the  mesosphere  (results  not  shown). 

In  addition  to  having  a  strong  impact  on  climato¬ 
logical  mean  quantities,  parametrized  OGWD  also  has 
a  strong  impact  on  variability,  particularly  in  the  strato¬ 
sphere  (Siskind  et  al.,  2007;  Richter  et  al.,  2010).  Figure  2 
illustrates  the  annual  cycle  of  interannual  standard  devi¬ 
ation  of  temperature  at  80N  and  30  hPa  calculated  from 
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(a)  WEAK 


(b)  STRONG 


(c)  STRONG-WEAK 


Figure  1.  Zonal-mean  zonal  wind  (contours)  and  OGWD  (shading,  in  10  6  Pa  m  *)  in  NH  winter  (DJF)  for  control  simulations  with  a) 
weak  drag  (G= 0.25)  and  b)  strong  drag  (G=1.0),  and  c)  its  difference.  (After  Sigmond  and  Scinocca  (2010)). 


daily  data.  The  ERA40  reanalysis  (red)  shows  large  vari¬ 
ability  during  northern  winter  peaking  around  the  mid¬ 
dle  of  February,  comparatively  smaller  variability  during 
summer  and  a  progressive  increase  during  fall  into  early 
winter.  Results  from  two  simulations  of  WACCM  forced 
with  observed  SSTs  are  also  shown:  In  the  first  simula¬ 
tion  (dash-dot),  a  value  is  used  for  the  efficiency  factor 
e  that  produces  the  best  tropospheric  simulation  (Garcia 
et  ah,  2007);  in  the  second  simulation  (solid  black),  e  is 
reduced  by  a  factor  of  1 .6-1.  Overall,  both  simulations  are 
deficient  compared  to  the  reanalysis  when  contrasting  the 
standard  deviations  during  winter.  This  deficiency  leads 
to  a  reduced  number  of  stratospheric  warmings  (Richter 
et  ah,  2010).  However,  the  simulation  with  the  smaller 
value  of  e  (solid)  is  much  closer  to  the  reanalysis  and 
produces  more  realistic  statistics  of  stratospheric  warm¬ 
ings  (not  shown).  The  improved  variability  in  this  second 
model  simulation  is  due  to  small  but  important  changes 
of  the  zonal  mean  zonal  wind  in  the  lower  stratosphere, 
which  favor  a  more  poleward  concentration  of  planetary 
wave  activity  in  the  winter-time  polar  vortex  (not  shown). 


80N  -  30  hPa 


Figure  2.  Annual  cycle  of  northern  hemisphere  interannual  tem¬ 
perature  variability  determined  from  daily  data  from  ERA40  (red) 
and  two  WACCM  simulations  with  differing  values  of  the  OGWD 
efficiency  factor  e. 

Given  the  above  sensitivity  of  the  current  climate  to 
these  parameters,  it  is  important  to  quantify  their  impact 
on  climate  change  simulations.  Sigmond  et  al.  (2007) 
showed  that  in  the  GFDL  AM2  model  the  500  hPa  cir¬ 
culation  response  at  high  latitudes  depends  strongly  on 
the  tunable  parameter  G.  In  another  study,  Sigmond  et  al. 
(2008)  showed  that  the  Northern  Annular  Mode  (NAM) 
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response  to  a  CO2  doubling  perturbation  is  more  depen¬ 
dent  on  G  than  on  the  model  lid  height  in  the  CCCma 
model.  Their  result  left  open  the  question  of  whether  the 
NAM  response  to  CO2  doubling  (and  associated  regional 
climate  responses)  could  depend  directly  on  the  poorly 
constrained  parametrized  OGWD  response.  Alternatively, 
the  NAM  response  may  instead  be  mainly  sensitive  to  the 
model  basic  state  winds,  which  the  OGWD  parametriza- 
tion  settings  help  to  shape.  Sigmond  and  Scinocca  (2010) 
addressed  this  directly  with  controlled  model  experiments 
that  separately  doubled  CO2  and  varied  the  OGWD.  They 
found  that  the  basic  state  winds  were  the  controlling  factor 
shaping  the  climate  response  to  C02  doubling,  whereas 
the  parametrized  OGWD  response  had  very  little  direct 
effect.  The  important  effect  of  the  OGWD  is  to  minimize 
the  biases  in  the  lower  stratospheric  basic  state  and  to  set 
the  refractive  properties  for  resolved  waves. 

OGWD  has  a  strong  and  robust  impact  on  the 
climate-change  induced  increase  in  the  equator-to-pole 
Brewer-Dobson  circulation  in  the  stratosphere  in  models 
(Brewer,  1949)  as  inferred  from  Butchart  et  al.  (2006) 
and  confirmed  subsequently  by  Li  et  al.  (2008),  McLan- 
dress  and  Shepherd  (2009),  and  Butchart  et  al.  (2010).  The 
mechanism,  which  is  quite  straightforward,  is  as  follows: 
Increased  greenhouse  gases  warm  the  troposphere  and 
cool  the  stratosphere.  This  increases  the  meridional  tem¬ 
perature  gradient  in  the  subtropical  middle  and  upper  tro¬ 
posphere,  which,  through  thermal  wind  balance,  increases 
the  strength  of  westerlies  on  the  upper  flank  of  the  sub¬ 
tropical  jet.  Since  a  large  fraction  of  the  parametrized 
orographic  waves  occur  over  the  Himalayas,  the  drag 
that  these  waves  produce  occurs  as  a  result  of  wave 
breaking  above  the  subtropical  jet  maximum.  As  climate 
change  increases  the  mean  wind  speeds  in  this  region, 
the  parametrized  waves  break  higher,  causing  increased 


wave  drag  in  the  lower  stratosphere.  (The  above  argu¬ 
ment  implicitly  assumes  that  OGWD  reacts  passively 
to  changes  in  the  zonal  wind  brought  about  by  climate 
change:  Sigmond  and  Scinocca  (2010)  demonstrated  that 
this  is,  in  fact,  the  case.)  Through  “downward  control”  this 
causes  increased  upwelling  on  the  equatorward  side  of  the 
OGWD  maximum. 

2.3  Issues  regarding  the  ozone  hole 

Because  parametrized  OGWD  depends  on  geographic 
variations  in  resolved  mesoscale  topography,  the  role  of 
OGWD  is  considerably  weaker  in  the  Southern  Hemi¬ 
sphere  than  in  the  Northern  Hemisphere.  Consequently, 
NGWD  has  played  a  larger  relative  role  in  controlling 
modelled  winds  and  temperatures  in  the  winter  and  spring 
in  the  Antarctic.  Garcia  and  Boville  (1994)  showed  that 
mesospheric  NGWD  can  significantly  impact  the  circula¬ 
tion  and  the  thermal  structure  in  the  stratosphere  through 
downward  control.  They  showed  that  the  effect  on  the 
temperature  and  vertical  residual  velocity  is  largest  in  the 
upper  stratosphere  and  decreases  at  lower  levels.  Although 
the  magnitudes  of  the  anomalies  associated  with  NGWD 
are  small  when  compared  to  the  mean  climatology,  chem¬ 
ical  processes  are  particularly  sensitive  to  slight  changes 
of  the  ambient  temperature  and  even  differences  of  a  few 
degrees  can  have  a  large  impact  on  the  activation  of  the 
heterogeneous  reactions.  In  this  respect,  the  timing  of 
those  changes  is  even  more  important.  In  fact,  Garcia 
and  Boville  (1994)  show  that  near  30  km  the  tempera¬ 
ture  difference  in  July  at  80S  can  be  as  large  as  10  K, 
when  comparing  simulations  with  and  without  NGWD; 
the  same  temperature  difference  in  September  is  smaller 
but  is  still  several  degrees  (see  their  Fig.  5).  Heteroge¬ 
neous  processes  critical  for  ozone  depletion  are  sensitive 
to  temperatures  in  the  lower  stratosphere  throughout  win¬ 
ter  and  early  spring  (Solomon,  1988). 
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Tuning  of  NGWD  parametrizations  with  interactive 
chemistry  presents  challenges  that  were  unforeseen  sev¬ 
eral  years  ago.  Without  interactive  chemistry,  tuning  was 
only  required  for  the  climatological  mean  temperatures 
and  zonal  mean  zonal  winds,  while  now  with  interactive 
chemisty,  tuning  for  the  seasonal  cycle  is  also  required. 
Figure  3  illustrates  the  difficulties  that  can  be  encoun¬ 
tered  with  interactive  chemistry.  Two  simulations  using 
WACCM  (Garcia  et  ah,  2007)  with  1995  boundary  condi¬ 
tions  for  chemical  composition  are  shown  for  the  south¬ 
ern  hemisphere.  The  two  simulations  differ  only  in  the 
amount  of  non-orographic  gravity  wave  momentum  flux  at 
the  source  level.  The  simulation  with  the  lower  value  (0.7 
mPa;  left  panels)  shows  a  more  realistic  southern  hemi¬ 
sphere  ozone  minimum.  By  contrast,  the  simulation  with  a 
larger  source  level  momentum  flux  (1.0  mPa;  right  panels) 
shows  a  much  weaker,  shorter  duration  ozone  hole. 

The  remarkable  difference  of  column  ozone  between 
the  two  simulations  of  Figure  3  is  related  to  seasonal 
temperature  variations  in  the  lower  stratosphere  which 
directly  impact  ozone  depletion  chemistry.  The  lower 
panels  of  Figure  3  illustrate  the  annual  march  of  lower 
stratospheric  temperature  at  61  hPa  over  the  south  polar 
cap  in  these  two  simulations.  With  larger  source  level 
momentum  flux  (lower  right)  a  more  rapid  increase  of 
temperature  occurs  during  southern  spring  compared  to 
consistently  lower  temperatures  between  September  and 
November  in  the  simulation  with  smaller  source  level 
momentum  flux  (lower  left).  Temperatures  in  the  left  panel 
remain  too  cold  in  December,  while  on  the  right,  the 
warming  in  the  lower  stratosphere  occurs  too  soon.  The 
October  monthly  mean  temperature  is  already  above  200 
K  in  the  simulation  with  larger  source  momentum  flux, 
sufficient  to  inhibit  the  formation  of  polar  stratospheric 
clouds  and  the  chemical  conversion  of  chlorine  into  active 


forms  that  deplete  ozone  in  early  spring.  The  two  settings 
of  source  level  momentum  flux  differ  by  only  a  small 
amount  relative  to  the  uncertainties  in  these  fluxes. 

3  Global  observations  of  gravity  wave  momentum 
fluxes 

The  above  model  issues  call  for  improved  constraints 
on  momentum  fluxes  in  the  lower  stratosphere.  The 
first  direct  observation  of  momentum  fluxes  (in  the 
mesosphere),  M  =  (Mx,  My)  =  p(u'w',  v'w'),  by  gravity 
waves  was  made  by  Vincent  and  Reid  (1983)  using  a 
split  beam  technique  with  an  FIF  radar.  This  technique 
has  since  been  used  in  the  stratosphere  at  VFIF  radar 
sites  in  Japan,  Peru,  the  UK,  and  India  beginning  in  the 
early  1990’s  (e.  g.  Sato,  1993;  Riggin  et  ah,  1997).  In 
recent  years  though,  observational  studies  using  several 
other  techniques  are  now  providing  global-scale  informa¬ 
tion  on  gravity  wave  momentum  fluxes  in  the  stratosphere. 
A  common  thread  in  these  observations  is  that  they  ana¬ 
lyze  fluctuations  about  some  mean  as  an  indicator  of  grav¬ 
ity  waves  and  then  estimate  the  associated  gravity  wave 
amplitudes,  and  use  gravity  wave  polarization  relations  to 
derive  the  momentum  fluxes  either  from  temperature  mea¬ 
surements  only  or  from  wind  and  temperature  measure¬ 
ments.  Other  than  radar  and  superpressure  balloon  studies, 
all  the  other  techniques  discussed  in  this  paper  use  spatial 
fluctuations  to  derive  gravity  wave  momentum  fluxes,  and 
all  techniques  have  inherent  horizontal  and  vertical  reso¬ 
lutions.  As  we  will  see,  these  spatial  resolution  limits  also 
imply  measurements  of  certain  gravity  wave  frequency 
ranges. 

3. 1  The  gravity  wave  dispersion  relation 

Gravity  waves  are  treated  here  and  in  the  analyses  cited  in 
the  paper  as  small-amplitude  perturbations  to  some  larger 
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Figure  3.  Top  panels:  Time-latitude  changes  in  southern  hemisphere  column  ozone  in  WACCM  simulations  with  slightly  lower  (left,  0.7 
mPa)  and  higher  (right,  1.0  mPa)  values  of  NGWD  source-level  momentum  flux.  Bottom  panels:  Time  series  of  polar  temperatures  at  the 

61  hPa  pressure  level  in  the  simulations. 


scale  horizontally  uniform  and  steady  background  state 
with  horizontal  wind  u,v,  temperature  T  and  density  p. 
We  refer  the  reader  to  Fritts  and  Alexander  (2003;  section 
2)  for  a  derivation  of  the  wave  equations  starting  from  the 
fundamental  fluid  equations,  which  can  be  used  to  derive 
the  dispersion  relation. 

The  dispersion  relation  relates  the  wave  frequency 
u>  to  other  wave  propagation  properties  such  as  the  hor¬ 
izontal  (k,  l )  and  vertical  (to)  wavenumbers.  An  important 
wave  frequency  that  emerges  is  the  intrinsic  frequency 
Gj  —  u  —  ku  —  Iv,  that  which  would  be  measured  in  the 
frame  of  reference  moving  with  the  background  wind.  The 
dispersion  relation  can  be  generally  written  (e.g.  Fritts  and 
Alexander  (2003);  their  eqn.  (23))  as 

N2(k2  +  l2)  +  /2[m2  +  (2H)~2] 
k2  +  l2  +  m2  +  (2H)~2 


where  N  is  the  Brunt- Vaisiila  frequency;  /  is  the  Cori¬ 
olis  parameter;  and  H  is  the  density  scale  height.  (We 
are  further  neglecting  compressibility  effects  that  lead  to 
acoustic-gravity  waves.)  With  slow  variations  in  height 
N(z),  H(z),  u(z),v(z),  ( 1 )  allows  variations  m(z)  and 
Gj{z).  Altitudes  where  Gj  — »  0  (or  c  =  u/k  =  U,  where  U 
is  the  wind  component  in  the  direction  of  wave  propaga¬ 
tion)  are  called  critical  levels.  Near  or  below  this  level 
the  wave  generally  becomes  unstable  and  will  dissipate. 
These  processes  underlie  the  concept  of  “critical  level  fil¬ 
tering”. 

Various  observation  techniques  have  different  inher¬ 
ent  resolutions,  and  the  method  of  separating  waves 
from  the  larger-scale  background  can  further  restrict  the 
range  of  wavelengths  that  are  resolved  in  an  analysis.  A 
restricted  range  of  wavelengths  also  implies  a  restricted 
range  of  frequencies  through  (1),  a  point  made  earlier 
by  Alexander  ( 1998).  We  will  describe  the  limitations  of 
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individual  techniques  in  sections  3. 3-3. 5  and  compare  the 
various  observable  ranges  of  wavelength  and  frequency  in 
section  3.7. 

3.2  Calculation  of  gravity  wave  momentum  fluxes 

The  method  of  calculating  gravity  wave  momentum  fluxes 
depends  on  which  variables  are  being  directly  observed. 
For  instance,  the  primary  satellite-observed  quantity  is 
temperature  derived  from  radiances.  In  this  case,  the 
gravity  wave  polarization  relations  can  be  used  to  derive 
the  implied  horizontal  and  vertical  wind  fluctuations, 
and  these,  in  turn,  can  be  used  to  derive  gravity  wave 
momentum  flux  using  the  following  formula  from  Ern  et 
al.  (2004): 


:*)’( 

fV 

(  k 

l  \ 

(MX,  My)  =  ~P  ( 

?) 
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\m 

-  (2) 
m  J 

where  (k,l)  =  (2tt/Xx, 

2n/Xy) 

is 

the 

horizontal 

wavenumber  vector,  m  =  27r/A2  is  the  vertical  wavenum¬ 
ber,  and  T  is  the  temperature  amplitude  of  the  wave, 
and  T  a  background  (or  larger-scale  mean)  temperature. 
Equation  (2)  uses  a  simplified  form  of  the  dispersion 
relation  appropriate  for  medium  frequency  waves, 
Cj  =  N\(k,l)/m\.  Radiosonde  techniques,  instead  keep 
terms  that  are  important  at  low  frequencies.  Appropriate 
assumptions  are  given  in  the  sections  that  follow. 

3.3  Radiosondes 

Standard  radiosonde  data  are  meant  to  supply  temper¬ 
ature,  wind,  and  humidity  data  for  weather  prediction 
models,  so  information  on  only  a  small  number  of 
mandatory  and  significant  levels  (e.g.  Eskridge  et  al., 
1995)  are  saved  and  sent  to  weather  prediction  centers. 
These  data  cannot  easily  be  used  to  obtain  informa¬ 
tion  on  gravity  waves.  Allen  and  Vincent  (1995)  illus¬ 
trated  how  gravity  wave  information  could  be  obtained 


using  all  measurements  in  the  original  radiosonde  profiles. 
This  motivated  an  effort  by  the  World  Climate  Research 
Program’s  (WCRPs)  Stratospheric  Processes  and  their 
Role  in  Climate  (SPARC)  project  to  save  these  high 
vertical-resolution  radiosonde  data  (Hamilton  and  Vin¬ 
cent,  1995).  For  instance,  there  are  presently  nine  years 
of  US  high  vertical-resolution  (nominally  30  m  verti¬ 
cal  resolution)  archived  on  the  SPARC  Data  Center  (see 
http  ://www.  spare  .sunysb.edu/html/hres  .html) . 

Radiosonde  data  contain  information  on  the  horizon¬ 
tal  and  vertical  wind  components  as  well  as  temperature, 
but  Geller  and  Gong  (2009)  illustrated  how  different  grav¬ 
ity  wave  frequencies  are  preferentially  observed  in  these 
different  variables.  For  this  reason,  gravity  wave  momen¬ 
tum  fluxes  are  derived  from  the  high  vertical-resolution 
radiosonde  data  using  the  wind  and  temperature  measure¬ 
ments  together  with  the  gravity  wave  polarization  rela¬ 
tions  and  the  assumption  that  wave  energy  flux  is  upward 
by  the  following  formulae, 

(Mx,  My)  =  -p^(u'f+90,v'f+90)5.(uj)  (3) 

where  <5_(w)  =  (1  —  f2 /Cj2)  and  T+90  is  determined  by 
Hilbert  transforming  T  by  90°  and  then  taking  the  reverse 
transform  (see  Vincent  et  al.,  1997).  Using  the  US  high 
vertical-resolution  radiosonde  data,  Gong  et  al.  (2008) 
show  time  series  of  the  vertical  fluxes  of  zonal  and 
meridional  momentum  by  gravity  waves  at  individual 
radiosonde  launch  sites.  These  results  averaged  in  longi¬ 
tude  over  all  sites  are  shown  in  Figure  4. 

The  zonal  momentum  fluxes  tend  to  be  positive  (east¬ 
ward)  at  low  latitudes,  negative  (westward)  at  middle  lat¬ 
itudes,  and  positive  during  the  summer  months  at  high 
latitudes.  This  pattern  is  broadly  consistent  with  what 
would  be  expected  from  critical  level  filtering  by  the 
background  winds:  There  are  primarily  westward  winds 
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Time  series  of  zonal  momentum  flux  (pa)*  1000  in  Lower  Stratosphere 
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Time  series  of  meridional  momentum  flux  (pa)*  1000  in  Lower  Stratosphere 
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Figure  4.  Time-latitude  plots  of  monthly  mean  eastward  (top) 
and  northward  (bottom)  momentum  fluxes  derived  from  US  high 
vertical-resolution  radiosondes  in  units  of  lO”3  Pa  for  the  lower 
stratosphere  (18-25  km).  The  9-yr  time  series  at  each  latitude  is 
an  average  of  sites  in  that  latitude  bin.  Locations  of  the  sondes  are 
shown  in  Gong  et  al.  (2008).  These  are  extensions  of  earlier  results 
by  Wang  (2003). 


at  the  low  latitudes  in  the  upper  troposphere/lowermost 
stratosphere,  and  the  westward  winds  extend  furthest  pole- 
ward  in  mid-summer.  At  middle  latitudes  the  eastward 
jet  stream  winds  dominate  the  filtering,  giving  rise  to 
the  negative  zonal  momentum  fluxes  there.  The  picture 
for  the  meridional  momentum  fluxes  is  less  coherent,  but 
generally,  the  momentum  fluxes  are  equatorward  at  low 
latitudes,  poleward  in  middle  latitudes,  and  equatorward 
again  at  high  latitudes,  possibly  a  signature  of  the  three- 
celled  meridional  circulation  of  the  troposphere.  As  indi¬ 
cated  earlier,  these  radiosonde-derived  momentum  fluxes 
are  for  gravity  waves  with  low  intrinsic  frequencies,  com¬ 
monly  Cj  ~  2/  —  4/,  and  horizontal  wavelengths  ~500- 
2000  km. 


3.3.1  Long-duration  balloons 

Superpressure  balloons  are  closed  and  made  of  an  almost 
inextensible  material,  so  that  they  maintain  a  fixed  con¬ 
stant  volume  during  their  flights  provided  the  internal  gas 
pressure  exceeds  that  of  the  atmosphere.  Hence  these  bal¬ 
loons  drift  on  constant-density  (isopycnic)  surfaces  in  the 
atmosphere,  and  are  capable  of  long-duration  flights  (up 
to  several  months)  that  cover  wide  geographical  areas 
(e.g.  The  TWERLE  Team,  1977;  Hertzog  et  ah,  2007). 
Observations  of  meteorological  variables  made  at  regular 
time  intervals  along  the  flight  permit  detection  of  dis¬ 
turbances  associated  with  gravity  waves.  A  unique  fea¬ 
ture  of  long-duration  balloon  observations  is  that  they  are 
collected  in  the  intrinsic  reference  frame:  Balloons  are 
advected  by  the  horizontal  wind,  and  therefore  behave 
as  quasi-Lagrangian  tracers  of  horizontal  motion  so  that 
wave  intrinsic  frequencies  are  directly  measured. 

Gravity-wave  momentum  flux  calculations  based 
on  long-duration  balloon  observations  were  pioneered 
by  Massman  (1981)  using  results  from  a  few  balloons 
launched  during  the  Tropical  Wind,  Energy  Conversion, 
and  Reference  Level  Experiment  (TWERLE)  in  1975- 
76.  Massman’s  results  were,  however,  confined  to  short 
periods  of  time  because  of  the  limited  data-transmission 
capabilities  available  at  that  time.  These  limitations  are 
no  longer  an  issue  because  of  advances  in  space-borne 
communications  and  positioning  systems.  Estimations  of 
gravity-wave  momentum  fluxes  representative  of  equato¬ 
rial  latitudes  were  derived  with  3  long-duration  balloons 
flying  in  the  lower  stratosphere  (around  60  hPa)  in  1998 
(Hertzog  and  Vial,  2001).  Recently,  global  momentum- 
flux  maps  were  obtained  at  the  same  level  over  the  poles  of 
both  hemispheres  with  the  dataset  collected  in  the  frame¬ 
work  of  the  Strateole/Vorcore  project  (Vincent  et  ah,  2007; 
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Hertzog  et  al.,  2008).  In  these  latter  studies,  GW  momen¬ 
tum  fluxes  are  computed  from  the  correlation  between 
the  horizontal  velocity  and  atmospheric  pressure  fluctu¬ 
ations.  The  pressure  disturbance  measured  by  the  bal¬ 
loon  is  linked  to  the  wave-induced  vertical  velocity  dis¬ 
turbance  (Boccara  et  al.,  2008).  It  was  also  assumed  that 
the  balloons  behave  as  perfect  isopycnic  tracers,  which 
is  a  realistic  approximation  for  the  waves  examined  in 
these  studies  but  may  be  not  applicable  for  shorter  period 
waves  (Nastrom,  1980)  that  will  be  resolved  in  future 
campaigns.  For  Strateole/Vorcore,  the  data  rate  was  lim¬ 
ited  to  (15  min)-1,  so  waves  with  intrinsic  periods  >1 
hr  were  resolved.  The  range  of  wavelengths  that  can  be 
analyzed  with  the  existing  data  will  be  further  described 
in  section  3.7  and  illustrated  graphically  in  Figure  8  in 
context  with  other  measurement  techniques.  The  hydro¬ 
static  approximation  to  the  observational  limit  of  these 
data  gives  k  <  m/12. 

Results  of  the  Vorcore  campaign  show  zonal-mean 
high-latitude  momentum  fluxes  ~  2  —  3  mPa,  but  larger 
average  values  near  30  mPa  in  a  10°  x  5°  region  in 
the  vicinity  of  the  Antarctic  Peninsula.  The  fluxes  are 
largely  westward,  but  significant  localized  meridional 
fluxes  also  occur  over  topography.  The  intermittency  in 
these  measured  fluxes  will  be  described  in  section  3.6. 

3.4  Nadir  scanning  satellites 

For  the  past  30  years,  nadir-viewing  passive  microwave 
and  infrared  remote  sensors  have  been  a  cornerstone 
of  operational  numerical  weather  prediction,  particularly 
thermal  02  and  CO2  channels  from  which  temperatures 
can  be  derived.  However,  until  recently,  two  properties 
of  these  measurements  made  them  incapable  of  resolv¬ 
ing  gravity  waves.  First,  the  nadir  view  yields  broad 
weighting  functions  in  the  vertical,  such  that  vertically 
short  gravity-wave  oscillations  are  averaged  out  in  these 


nadir  radiances.  Second,  early  instruments  had  horizontal 
resolutions  of  several  hundred  kilometers,  making  them 
unable  to  resolve  typical  gravity  wave  horizontal  wave¬ 
length  scales  of  ~5-500  km.  For  example,  the  cross-track 
scanning  Microwave  Sounding  Unit  (MSU),  which  flew 
on  the  NOAA-6  through  NOAA-14  satellites,  had  hori¬ 
zontal  measurement  footprints  in  excess  of  100  km. 

Starting  with  the  NOAA-15  satellite,  launched  in 
1998,  MSU  was  superseded  by  the  Advanced  MSU  Unit- 
A  (AMSU-A),  whose  footprint  diameters  were  reduced 
to  --40  km  at  nadir  and  whose  radiometric  precision 
was  improved  to  ~0.2-0.8  K  in  the  middle  atmosphere. 
Wu  (2004)  isolated  fluctuations  in  along-track  AMSU-A 
radiances  that  revealed  enhancements  well  above  nomi¬ 
nal  noise  floors  near  the  vortex  edge  and  over  orography 
that  correlated  with  similar  gravity-wave-induced  radi¬ 
ance  enhancements  observed  by  the  Upper  Atmosphere 
Research  Satellite-Microwave  Limb  Sounder  (UARS- 
MLS).  The  launch  of  NASA’s  Aqua  satellite  in  2002 
saw  AMSU-A  integrated  with  the  Atmospheric  Infrared 
Sounder  (AIRS)  to  form  a  high-resolution  coupled  tem¬ 
perature  sounding  system  (Parkinson,  2003).  In  addition 
to  many  more  thermal  channels  for  analysis,  the  horizon¬ 
tal  footprint  for  the  AIRS  measurements  are  three  times 
smaller  than  the  corresponding  AMSU-A  measurements, 
while  channel  noise  levels  are  comparable.  This  improved 
horizontal  resolution  makes  AIRS  sensitive  to  gravity 
waves  of  shorter  horizontal  wavelength  than  AMSU-A. 

AMSU-A  and  AIRS  are  not  purely  nadir  viewers: 
The  instruments  scan  cross-track  symmetrically  about 
nadir.  When  coupled  with  the  satellite  motion,  a  two- 
dimensional  radiance  image  with  a  cross-track  swath 
width  of  ~2000  km  is  assembled.  Wu  and  Zhang  (2004) 
presented  initial  evidence  that  horizontal  wave  structure 
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was  captured  in  AMSU  imagery.  To  investigate  this,  Eck- 
ermann  and  Wu  (2006),  using  a  simple  in-orbit  forward 
model  of  lower  stratospheric  (channel  9)  radiance  acquisi¬ 
tion,  showed  that  two-dimensional  horizontal  phase  struc¬ 
ture  of  large- amplitude  long-wavelength  gravity  waves 
could  be  resolved  (e.g.,  horizontal  wavelength  and  prop¬ 
agation  direction).  Eckermann  et  al.  (2006)  validated  the 
model  predictions  for  a  large-amplitude  mountain  wave 
over  Scandinavia,  reproducing  the  observed  AMSU-A 
radiance  perturbations  via  forward  modelling  of  high- 
resolution  numerical  weather  prediction  model  simula¬ 
tions  of  the  event.  Eckermann  et  al.  (2007)  extended  the 
analysis  of  that  wave  to  all  AMSU-A  stratospheric  chan¬ 
nels  9-14,  which  revealed  the  wave  propagating  right 
through  the  stratosphere  and  evolving  in  phase  structure 
with  height  due  to  directional  wind  shear. 

Alexander  and  Barnet  (2007)  showed  large- 
amplitude  mountain  wave  oscillations  imaged  in  AIRS 
swath  radiances  from  an  upper- strata  spheric  15  ji  m 
CO2  channel.  Subsequent  studies  of  radiances  from 
various  stratospheric  AIRS  channels  have  revealed  large- 
amplitude  stratospheric  mountain  waves  over  various 
mountainous  regions  (e.g.,  Alexander  and  Teitelbaum, 
2007;  Eckermann  et  al.,  2007;  Limpasuvan  et  al.,  2007; 
Eckermann  et  al.,  2009;  Alexander  et  al.,  2009)  as  well 
as  stratospheric  gravity  waves  generated  by  convection 
(Grimsdell  et  al.,  2009).  Gravity  waves  imaged  in  AIRS 
4.3//m  thermal  channels  validate  earlier  imaging  of 
apparent  upper-stratospheric  gravity  wave  structures  in 
these  band  radiances  with  the  SPIRIT-3  radiometer  (e.g., 
Picard  et  al.,  1998;  Dewan  et  al.,  1998). 

To  date,  studies  have  focused  primarily  on  radiances, 


rather  than  temperature  retrievals,  due  to  the  closer  con¬ 
nection  between  the  measured  perturbation  and  the  prop¬ 
erties  of  the  gravity  wave  that  produced  it  (e.g.,  Eck¬ 
ermann  and  Wu,  2006).  Nonetheless,  an  accurate  wave- 
induced  temperature  amplitude  is  required  to  deduce  more 
derived  wave  properties,  such  as  momentum  flux.  One 
approach  is  to  derive  first  the  wave’s  vertical  and  hori¬ 
zontal  wavelengths  from  the  radiance  data,  then  use  those 
wavelengths  to  quantify  the  vertical  and  horizontal  smear¬ 
ing  of  the  wave  due  to  the  vertical  weighting  functions  and 
horizontal  footprints,  respectively,  and  then  use  that  infor¬ 
mation  to  correct  the  radiance  perturbation  amplitude  for 
this  smearing  to  derive  the  wave’s  true  temperature  ampli¬ 
tude.  This  process  works  well  in  case  studies  (e.g.,  Ecker¬ 
mann  et  al.,  2006;  Alexander  and  Teitelbaum,  2007;  Eck¬ 
ermann  et  al.,  2009;  Alexander  et  al.,  2009),  but  must  be 
done  individually  and  carefully  for  each  event  in  question. 
At  lower  altitudes  where  AIRS  radiances  become  affected 
by  clouds,  this  fails.  Here,  a  detailed  AIRS/AMSU-A  tem¬ 
perature  retrieval  algorithm  is  needed  in  which  AMSU-A 
radiances  are  used  to  ’’cloud-clear”  AIRS  radiances  prior 
to  retrieving  temperatures  (Susskind  et  al.,  2006).  These 
AIRS/AMSU-A  temperature  retrievals  sometimes  retain 
gravity  wave  structures  seen  in  the  original  radiances 
(e.g.,  Alexander  and  Teitelbaum,  2007;  Eckermann  et  al., 
2007)  but  at  the  expense  of  desampling  to  the  coarser 
AMSU-A  footprint  resolution.  Whether  gravity  wave  ver¬ 
tical  structure  is  retrieved  optimally  by  these  algorithms 
is  not  currently  clear.  Recently,  Hoffmann  and  Alexan¬ 
der  (2009)  developed  an  AIRS-only  temperature  retrieval. 
While  this  retrieval  preserves  the  fine  horizontal  resolu¬ 
tion  of  the  original  AIRS  data,  it  has  higher  noise  and 
does  not  extend  to  lower  altitudes  like  the  AIRS/AMSU-A 
temperature  retrievals.  Nonetheless,  for  short  deep  gravity 
waves  that  are  believed  to  carry  much  of  the  gravity- wave 
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momentum  flux,  the  Hoffmann- Alexander  AIRS  tempera¬ 
ture  retrievals  look  to  be  an  exciting  new  product  to  better 
quantify  the  gravity  wave  momentum  budget  globally  in 
the  middle  atmosphere. 

Momentum  fluxes  have  been  derived  in  the  win¬ 
ter  stratosphere  for  mountain  wave  events  observed  over 
Scandinavia,  the  Antarctic  Peninsula,  and  South  Georgia 
Island  with  AMSU  and  AIRS  using  equation  (2).  Peak 
values  in  the  Scandinavian  event  were  300  mPa  (Eck- 
ermann  et  ah,  2007)  and  for  the  Antarctic  event  140 
mPa  (Alexander  and  Teitelbaum,  2007).  These  were  wave 
events  with  horizontal  wavelengths  of  400  and  300  km 
respectively.  Events  over  S.  Georgia  described  in  Alexan¬ 
der  et  al.  (2009)  included  a  range  of  horizontal  wave¬ 
lengths  from  ~50-400  km.  Averaged  over  a  2.5°  x  2.5° 
area,  the  S.  Georgia  momentum  fluxes  ranged  60-200 
mPa,  while  local  values  were  >1000  mPa.  Despite  these 
large  fluxes,  the  island  is  too  small  to  be  identified  as  a 
land  point  in  global  climate  models,  and  the  parameter¬ 
ized  OGWD  is  zero.  Fluxes  from  islands  in  the  southern 
oceans  may  be  an  important  missing  source  of  wave  drag 
for  the  southern  hemisphere  circulation. 

3.5  Limb  scanning  satellites 

Infrared  limb  sounding  was  the  first  technique  used  to 
detect  GWs  from  space  (Fetzer  and  Gille,  1994).  In 
this  geometry  the  instrument  views  slightly  downward 
towards  the  Earth's  horizon,  through  the  atmosphere  and 
into  cold  space  (Preusse  et  ah,  2002).  The  instrument 
measures  the  infrared  radiation  from  thermal  emissions 
along  the  line  of  sight  (LOS).  Most  of  the  radiance  stems 
from  the  lowest  part  of  the  LOS  close  to  the  tangent 
height,  where  the  LOS  is  almost  parallel  to  the  Earth's 
surface.  Therefore  infrared  limb  sounding  is  sensitive 
to  GWs  with  vertical  wavelengths  longer  than  2-5  km 
(depending  on  instrument),  and  primarily  sensitive  to 


GWs  of  horizontal  wavelengths  longer  than  100-200  km 
(Preusse  et  ah,  2002).  Altitude  profiles  are  taken  by 
varying  the  viewing  angle  and  thus  the  tangent  height. 
Typically,  the  registration  of  one  altitude  profile  by  an 
infrared  limb  sounder  takes  between  10s  and  one  minute, 
corresponding  to  80  -  500  km  distance  separation  of 
adjacent  profiles  along  the  orbit  track. 

Infrared  limb  sounders  have  been  flown  on  high- 
inclination  low  Earth  orbits,  resulting  in  typical  orbit  dura¬ 
tions  of  90-100  minutes  and  ~  1 4  orbits  per  day.  Measure¬ 
ments  during  day  and  night  allow  for  precise  estimates 
of  the  zonal  mean  background  and  planetary  scale  waves 
with  global  wavenumbers  up  to  ~7  (Salby,  1982).  The 
background  atmosphere  (estimated  as  zonal  mean  plus 
global-scale  waves)  is  subtracted  from  the  data  and  the 
residual  temperature  fluctuations  are  interpreted  as  gravity 
waves.  Many  case  studies  and  global  investigations  show 
that  GWs  can  be  reliably  inferred  in  this  way. 

For  typical  satellite  resolution,  the  momentum  flux 
can  be  approximated  with  equation  (2).  Applying  wave 
analyses  to  the  vertical  profiles  of  temperature  residu¬ 
als  we  can  infer  the  vertical  wavelengths  of  the  most 
prominent  wave  signatures  and  the  corresponding  tem¬ 
perature  amplitudes.  The  horizontal  wavelength,  however, 
is  undersampled  (Ern  et  al.,  2004).  In  addition,  GWs 
in  the  lower  stratosphere  tend  to  occur  in  packets  with 
limited  horizontal  extent.  A  good  example  is  the  moun¬ 
tain  wave  event  above  the  Andes  reported  by  Eckermann 
and  Preusse  (1999).  Three  consecutive  profiles  above  the 
Andes  showed  distinct  wave  signatures  with  the  proper¬ 
ties  of  mountain  waves.  The  profiles  measured  before  and 
afterwards  over  the  ocean  do  not  show  pronounced  wave 
signatures.  Because  of  this  limited  extent,  an  innovative 
way  of  determining  the  horizontal  wavelength  is  required: 
A  common  dominant  wave  signal  is  identified  in  a  pair 
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The  CRISTA  instrument  was  deployed  on  two  one- 
week  long  Space  Shuttle  missions  in  November  1994 
and  August  1997.  Vertical  profiles  were  investigated  by  a 
combination  of  Maximum  Entropy  Method  and  sinusoidal 
fit.  The  momentum  flux  was  calculated  from  the  average 
vertical  wavelength  and  temperature  of  the  leading  wave 
component  of  adjacent  profiles  (for  details  see  Ern  et 
al.,  2004).  The  momentum  flux  estimates  were  corrected 
for  the  instrumental  sensitivity  and  aliasing.  Results  are 
reproduced  in  Figure  5  a.  Correction  factors  have  been 
applied  to  these  “best  estimate”  values  to  approximately 
account  for  the  low  bias  associated  with  unknown  wave 
propagation  directions  (Ern  et  al.,  2004).  The  results  show 
very  similar'  longitudinal  variations  as  with  momentum 
fluxes  deduced  from  superpressure  balloons  (Hertzog  et 
al.,  2008). 


Figure  5.  Global  maps  of  gravity  wave  momentum  flux  at  20- 
30  km  altitude  derived  from  CRISTA  (top)  and  HIRDLS  (bottom) 
high  resolution  satellite  temperature  observations.  CRISTA  data  for 
this  map  were  obtained  Aug.  8-16,  1997,  and  HIRDLS  data  were 
obtained  Aug.  1-31,  2006.  After  Ern  et  al.  (2004)  and  Alexander  et 
al.  (2008) 


The  High  Resolution  Dynamics  Limb  Sounder 
(HIRDLS),  on  the  EOS-Aura  satellite,  provided  high 
vertical  resolution  temperature  measurements  from  Jan. 
2005  to  Mar'.  2008.  Vertical  profiles  were  analyzed  by 
means  of  the  S-transform  and  the  cross-spectrum  ampli- 


of  adjacent  vertical  profiles,  and  horizontal  wavenumber 
along  the  orbit  track  is  inferred  from  the  difference  of  the 
phase  ^  between  the  two  profiles  i  and  j  with  spacing 
A Xij  (Ern  et  al.,  2004): 


_  d(j){xh)  _  Afaj 
h  dxh  A  Xid 


(4) 


tude  between  profile-pairs  was  used  to  identify  the  dom¬ 
inant  wave  mode  as  a  function  of  height  in  the  profiles 
(Alexander  et  al.,  2008).  The  horizontal  wavelength  was 
then  determined  using  the  method  in  (4)  and  shown  in 
Figure  5b.  No  “best  estimate”  corrections  were  applied 
to  the  HIRDLS  data  so  they  represent  lower  limits  of  the 
momentum  flux. 


Here  kh  is  the  component  of  the  total  horizontal  wavenum¬ 
ber  resolved  along  the  line  joining  the  two  profiles.  Using 
this  method  global  momentum  flux  estimates  have  been 
deduced  from  data  obtained  by  two  orbiting  instruments: 
The  Cryogenic  Infrared  Spectrometers  and  Telescopes 
for  the  Atmosphere  (CRISTA)  and  the  High  Resolution 
Dynamics  Limb  Sounder  (HIRDLS). 


Both  datasets  agree  well  in  their  general  structure. 
Largest  values  are  found  in  the  southern  polar  vortex.  A 
second  maximum  is  found  above  subtropical  convection 
in  the  summer  hemisphere.  In  the  tropics,  long  horizontal 
wavelength  GWs  are  dominant,  resulting  in  relatively  low 
momentum  flux  values.  At  mid  and  high  latitudes  of  the 
summer  hemisphere,  weak  stratospheric  winds  and  a  wind 
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reversal  between  troposphere  and  stratosphere  strongly 
reduce  the  GW  momentum  flux. 

In  absolute  values  the  two  results  are  different  by  a 
factor  of  ~5.  CRISTA  values  are  larger,  due  partly  to  the 
sensitivity/aliasing  corrections  applied  in  Ern  et  al.  (2004), 
partly  to  the  longer  vertical  wavelengths  included  in  the 
CRISTA  analysis,  and  as  discussed  in  the  next  section, 
partly  to  intermittency  and  the  different  averaging  times 
( 1  wk  vs  1  mo).  By  further  improving  the  wave  analysis 
method,  calibrating  the  inferred  values  to  other  methods, 
and  inferring  improved  corrections,  new  best  estimates 
can  be  obtained  in  the  future. 

3.6  Quantifying  wave  intermittency 

Gravity  waves  are  intermittent  in  nature.  Observations 
show  they  often  occur  in  isolated  large  amplitude 
wave  packets.  For  example,  the  August-mean  map  from 
HIRDLS  observations  in  Fig  5b  shows  maximum  fluxes 
of  6  mPa,  however  individual  measurements  exceed  100 
mPa.  These  large  amplitude  events  move  from  location  to 
location  on  different  days,  resulting  in  the  smaller  monthly 
average.  Such  behaviour  is  furthermore  exemplified  by 
gravity-wave  momentum-flux  probability  density  func¬ 
tions  (pdf).  Figure  6  for  instance  displays  momentum-flux 
pdfs  obtained  with  observations  collected  by  HIRDLS 
between  50S  and  64S  and  by  long-duration  balloons  south 
of  50S  during  Vorcore.  The  sporadic  character  of  gravity- 
wave  activity  is  linked  to  the  broad  tail  consistently 
exhibited  by  both  pdfs,  which  corresponds  to  relatively 
rare  events  with  large  momentum- flux  values.  As  shown 
on  this  figure,  these  large  values  can  exceed  the  mean 
momentum  flux  by  more  than  one  order  of  magnitude.  The 
differences  between  the  pdfs  likely  arise  because  of  dif¬ 
ferences  in  space  and  time  sampling  (August  for  HIRDLS 
and  no  observation  above  Antarctica,  September-October 
for  Vorcore  balloons  with  observations  above  Antarctica), 


Momentum— flux  pdf 


Figure  6.  Gravity-wave  momentum-flux  probability  density  func¬ 
tions  obtained  with  (bold  gray)  HIRDLS  and  (light  black)  long- 
duration  balloon  observations  (see  text  for  more  details).  The  mean 
value  for  each  pdf  is  indicated  by  the  arrow  at  the  bottom  of  the 
figure. 

and  also  differences  in  detection  of  small  amplitude  grav¬ 
ity  waves  in  both  datasets,  which  impact  the  lowest  bin  of 
the  distribution. 

Average  values  of  wave  amplitude  can  furthermore 
depend  on  the  area  and  time  over  which  the  average  is 
taken.  Ideally,  when  compared  with  parametrized  gravity 
waves,  observational  intermittency  like  this  would  thus 
need  to  be  quantified  on  the  same  grid  as  the  global 
model  using  the  parametrization.  An  additional  problem 
is  that  waves  from  orographic  and  nonorographic  sources 
cannot  generally  be  clearly  separated  for  computation  of 
the  necessary  OGWD  and  NGWD  efficiency  factors. 

Hertzog  et  al.  (2008)  nevertheless  presented  geo¬ 
graphical  information  on  gravity-wave  intermittency  in 
the  Vorcore  super-pressure  balloon  observations.  Intermit¬ 
tency  in  this  work  was  quantified  via  e\  =  (1  +  ct2//t2)-1 
where  //  is  the  mean  momentum  flux  in  each  geographical 
grid  box,  and  a2  the  variance  (cf.  Figure  7).  The  result  is 
labeled  “probability  of  observation”,  where  a  value  of  1 
represents  no  intermittency,  and  a  small  value  represents  a 
high  degree  of  intermittency. 

What  stands  out  in  Fig.  7  is  that  the  most  intermittent 
sources,  associated  with  the  lowest  probability  of  obser¬ 
vation,  are  found  above  mountainous  areas,  such  as  the 
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pa<u'//W'>  intermittency:  Bernoulli 


Figure  7.  Wave  intermittency  quantified  as  “probability  of  observa¬ 
tion”  computed  from  super-pressure  balloon  data  obtained  during 
the  Vorcore  campaign.  (Hertzog  et  al.  (2008)) 


dissipation  with  height  is  included  explicitly.  For  OGWD 
parametrizations,  intermittency  in  sources  is  also  included 
explicitly  through  the  dependence  of  the  wave  amplitude 
on  the  local  topography  and  the  surface  wind  and  sta¬ 
bility.  Some  source  intermittency  is  also  explicit  in  the 
newer  source  parametrizations  for  wave  emission  from 
convection  and  fronts  (e.g.  Richter  et  ah,  2010).  NGWD 
parametrizations  that  are  not  physically  linked  to  specific 
generation  mechanisms  lack  this  part  of  the  intermittency. 
Both  OGWD  and  NGWD  parametrizations  also  include 
dimensionless  scaling  factors,  such  as  G  and  e  (Section 
2.2)  that  are  sometimes  called  intermittency  factors.  Con¬ 
straints  on  these  parametrization  factors  will  require  care¬ 
ful  case  study  comparisons. 


Antarctic  Peninsula,  the  Ellsworth  Range,  and  the  Antarc¬ 
tic  coast  in  several  places,  including  Adelie  Land.  In  con¬ 
trast,  the  probability  of  observing  gravity  waves  is  largest, 
and  the  intermittency  smallest,  above  the  oceans  and  to  a 
lesser  extent  above  the  Antarctic  Plateau.  Hertzog  et  al. 
(2008)  also  discussed  alternate  methods  of  defining  inter¬ 
mittency  such  as  ratio  of  the  50%  percentile  (the  median) 
to  the  90%  percentile.  The  values  so  obtained  (not  shown) 
differ  from  Figure  7  in  magnitude,  but  the  patterns  remain 
very  similar. 

The  intermittency  seen  in  observations  arises  from 
two  causes:  Sources  of  the  waves  are  intermittent,  and 
waves  may  be  filtered  as  they  propagate  through  the  inter¬ 
vening  atmosphere  between  their  source  and  the  obser¬ 
vation  point.  There  are  also  two  types  of  filtering:  One 
where  the  wave  energy  (and  momentum  flux)  disappears 
because  of  dissipation,  and  the  other  where  the  wave  may 
be  refracted  by  changes  in  wind  and  stability  to  a  vertical 
wavelength  lying  outside  the  observable  range. 

Gravity  wave  parametrizations  also  represent  inter¬ 
mittency  in  several  ways.  The  intermittency  due  to  wave 


3.7  Mapping  of  observations  into  wavenumber- 
frequency  space 

Internal  atmospheric  gravity  waves  have  intrinsic  frequen¬ 
cies  between  the  buoyancy  frequency  and  the  Coriolis  par¬ 
ameter,  horizontal  wavelengths  ranging  from  a  few  kilo¬ 
meters  to  a  few  thousand  kilometers  and  vertical  wave¬ 
lengths  can  theoretically  vary  from  0  to  oo.  No  single 
instrument  can  cover  all  these  scales.  However,  to  inter¬ 
compare  fluxes  derived  from  different  methods  or  models 
in  order  to  quantify  the  contributions  of  parts  of  the  spec¬ 
trum  to  the  whole,  it  is  necessary  to  quantify  the  observa¬ 
tional  range  of  each  individual  instrument.  If  this  observa¬ 
tional  filter  is  known,  it  can  then  be  applied  to  model  data 
and  meaningful  comparisons  can  be  performed.  Studies 
have  shown  that  the  observational  filter  not  only  changes 
the  absolute  value  but  also  the  relative  global  distributions 
(Alexander,  1998;  Jiang  et  al,  2004;  Ern  et  al.  2004;  2006; 
Choi  et  al.  2009). 

It  is  important  to  keep  in  mind  that  the  observational 
filter  is  not  dependent  on  the  measurement  instrument 
alone.  The  first  step  in  analyzing  a  dataset  for  GWs  is 
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to  isolate  the  GWs  from  other  atmospheric  signals,  such 
as  planetary  scale  waves.  This  filtering  must  be  taken 
into  account  when  calculating  the  observational  filter.  For 
instance,  radiosondes  are  sensitive  to  almost  all  scales  of 
GWs.  Flowever,  GWs  are  determined  from  single  profiles 
by  a  high-pass  filter  that  rejects  all  vertical  wavelengths 
longer  than  ~7  km.  For  infrared  limb  sounding,  on  the 
other  hand,  the  global  background  can  be  estimated  and 
the  visibility  filter  is  defined  by  the  radiative  transfer  and 
the  vertical  field  of  view  of  the  instrument  (Preusse  et  al., 
2002). 

Depending  on  the  type  of  measurement,  the  visibility 
filter  is  best  defined  either  in  horizontal  wavelength  or 
in  frequency  space.  Figure  8  shows  both  representations. 
The  visibility  filters  have  been  converted  based  on  the 
GW  dispersion  relation  neglecting  the  Coriolis  force,  but 
including  corrections  for  nonhydrostatic  effects: 


period  [hr] 

5  2  1  0.5 


W  k2  +  m2  +  1/AH2  ’ 

The  conversion  used  typical  stratospheric  values  of 
N  =0.02  s-1  for  the  buoyancy  frequency  and  //  =6.5  km 
for  the  scale  height. 

The  lines  in  Figure  8  give  typical  visibility  limits. 
Depending  on  the  wave  amplitude  and  structure  of  the 
background  atmosphere  these  sensitivity  limits  may  shift 
slightly.  It  should  also  be  noted  that,  in  particular,  for 
satellite  instruments  the  sensitivity  inside  the  given  visi¬ 
bility  limits  is  a  value  substantially  smaller  than  1 .  Dashed 
lines  in  Figure  8  indicate  waves  seen  only  at  favorable 
viewing  geometries.  A  detailed  discussion  of  the  visibility 
limits  of  various  satellite  instr  uments  is  given  by  Preusse 
et  al.  (2008).  The  shown  visibility  limits  are  based  on 
the  following  studies:  Infrared  limb  sounding  (Preusse  et 
al.,  2002;  adapted  with  HIRDLS  minimum  vertical  wave¬ 
length  limits),  microwave  sub-limb  (McLandress  et  al., 


Figure  8.  Typical  visibility  limits  as  functions  of  horizontal  and 
vertical  wavenumber  (top)  and  frequency/vertical  wavenumber  (bot¬ 
tom)  for  various  satellite  and  balloon  measurement  techniques. 
Shaded  regions  are  not  visible  to  any  of  the  techniques.  (After 
Preusse  et  al.,  2008.) 

2000;  adapted  for  a  2.5  km  vertical  field  of  view),  nadir 
viewing  satellites  (Eckermann  and  Wu,  2006;  adapted 
with  AIRS  minimum  horizontal  wavelength  limits),  super¬ 
pressure  balloons  (Hertzog  et  al.,  2008),  and  radiosondes 
(Wang  et  al.,  2005). 

Together,  this  set  of  measurements  provides  very 
good  spectral  coverage  of  the  momentum  flux  over  the  rel¬ 
evant  region  of  the  spectrum  as  shown  in  Fig.  8.  Gaps  in 
coverage  remain  for  waves  with  very  high  intrinsic  fre¬ 
quency,  although  these  may  be  partly  covered  by  future 
satellite  (Preusse  et  al.,  2008)  and  superpressure  balloon 
measurements.  A  collaborative  effort  among  researchers 
with  expertise  in  the  existing  measurements  could  lead  to 
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a  combined  set  of  accurate  constraints  for  parameteriza¬ 
tion  s. 

4  High  resolution  global  models 

Current  global  models  used  for  weather  prediction  and 
climate  projection  must  include  gravity  wave  parameter- 
izations  because  gravity  waves  are  usually  sub-grid  scale 
phenomena  or  under-resolved  in  both  the  horizontal  and 
vertical  in  such  models.  However,  recent  developments  in 
supercomputing  allow  us  to  simulate  gravity  waves  explic¬ 
itly  in  high-resolution  global  climate  models  (e.g.  Hamil¬ 
ton  et  ah,  1999;  Watanabe  et  ah,  2008).  Although  the  full 
spectral  range  of  gravity  waves  is  not  completely  resolved 
even  in  such  high-resolution  models,  recent  work  indi¬ 
cates  that  large-scale  wind  and  temperature  fields  can  be 
realistically  simulated  without  gravity  wave  parameteriza- 
tions  (Watanabe  et  ah,  2008).  This  means  that  the  momen¬ 
tum  balance  is  internally  retained  between  the  mean  fields 
and  disturbances  including  resolved  gravity  waves,  and 
suggests  the  behavior  of  the  resolved  gravity  waves  is 
also  realistic  in  terms  of  propagation  and  dissipation  in 
the  model  atmosphere.  Using  such  high-resolution  global 
models  the  effects  of  gravity  waves  can  be  quantitatively 
evaluated  as  functions  of  time  and  3D  space.  The  realism 
of  the  modelled  gravity  waves  can  be  confirmed  through 
comparison  with  observations  by  radars,  lidars,  radioson¬ 
des  and  satellites  (Sato  et  ah,  1999;  Watanabe  et  ah,  2008; 
Kawatani  et  ah,  2009). 

Of  course,  some  shortcomings  of  these  gravity-wave 
resolving  models  have  to  be  considered.  The  dynami¬ 
cal  characteristics  of  gravity  waves  depend  on  the  model 
equations,  e.g.,  hydrostatic  balance  is  assumed.  The  gen¬ 
eration  of  gravity  waves  from  convection  will  depend  on 
the  cumulus  parametrization  scheme.  Also,  previous  mod¬ 
elling  and  observational  work  shows  that  much  of  the 


momentum  flux  spectrum  remains  unresolved  in  such  high 
resolution  models  (Hamilton,  1996;  Hertzog  et  ah,  2001; 
Fritts  et  ah  2006).  In  this  section,  some  recent  results 
from  gravity-wave  resolving  global  model  simulations  are 
reviewed. 

4. 1  Description  of  gravity-wave  resolving  GCMs 

In  order  to  quantify  the  roles  of  gravity  waves  in 
the  global  momentum  budget,  gravity-wave  resolving 
(GWR)-GCMs  require:  1)  use  of  neither  gravity  wave 
parametrizations  nor  data  assimilation  of  observed  fields 
that  may  contaminate  explicitly-resolved  wave  motions, 
2)  no  extra  damping  of  wave  motions  such  as  that  typically 
used  in  numerical  weather  prediction  models,  3)  a  model 
top  near  the  0.01  hPa  level  or  higher  in  order  to  include 
mesospheric  GWD,  and  4)  a  spatial  resolution  sufficiently 
high  to  achieve  realistic  background  wind  and  thermal 
structures,  as  well  as  realistic  spatio-temporal  spectra  of 
resolved  gravity  waves.  Indeed,  GWR-GCMs  satisfying 
one  or  two  of  these  conditions  have  been  widely  used 
and  have  made  significant  contributions  to  gravity-wave 
statistics,  generation,  and  propagation  issues  (Miyahara 
et  ah,  1986;  OSullivan  and  Dunkerton,  1995;  Sato  et  ah, 
1999;  Sato,  2000;  Kawatani  et  ah,  2003,  2004,  2005;  Sato 
and  Yoshiki,  2008;  Kawatani  et  ah,  2009).  To  date,  the 
Geophysical  Fluid  Dynamics  Laboratory  (GFDL)  SKYHI 
GCM  (e.g.,  Hamilton  et  al,  1999;  Koshyk  et  ah,  1999) 
and  the  JAGUAR  GCM  (Watanabe  et  ah  2008;  Watanabe 
and  Miyahara,  2009)  have  a  model  top  near  0.01  hPa  or 
higher,  and  likely  satisfy  all  four  requirements.  Accord¬ 
ing  to  studies  with  these  models,  a  minimum  resolved 
horizontal  wavelength  (A hmin)  shorter  than  200  km  and 
a  vertical  layer  thickness  less  than  1  km  are  minimal 
requirements  in  order  to  achieve  realistic  simulation  of 
the  mean  winds  and  temperatures  throughout  the  middle 
atmosphere.  The  GWR-GCM  by  Watanabe  et  ah  (2008) 
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has  a  horizontal  spectral  resolution  of  T213  (A hmin  ~188 
km)  and  an  extremely  high  vertical  resolution  (~300  m)  in 
order  to  simulate  realistic  3D  propagation  of  gravity  waves 
in  sheared  mean  flows,  as  well  as  wave-mean  flow  inter¬ 
actions  in  the  vicinity  of  critical  levels.  Results  from  this 
model  are  briefly  reviewed  in  the  rest  of  section  4. 

4.2  Momentum  budgets 

Figure  9  shows  meridional  cross  sections  for  the  zonal 
mean  zonal  winds,  E-P  flux  vectors  and  E-P  flux  diver¬ 
gence  in  July,  which  are  simulated  in  the  GWR-GCM 
by  Watanabe  et  al.  (2008).  The  meridional  structures  of 
the  southern  hemisphere  polar  vortex  and  the  northern 
hemisphere  summertime  easterlies  in  the  stratosphere  and 
mesosphere  are  well  simulated.  The  wave  forcing  due  to 
simulated  gravity  waves  with  horizontal  wavelengths  of 
1 88-950  km  (GW  group)  dominates  the  total  wave  forcing 
in  the  mesosphere  of  both  hemispheres,  which  strongly 
decelerates  the  upper  part  of  the  polar  vortex  and  the  sum¬ 
mertime  easterlies  to  maintain  their  characteristic  merid¬ 
ional  structures.  The  meridional  distribution  of  the  E-P 
flux  in  Fig.  9  indicates  that  westward  propagating  grav¬ 
ity  waves  originating  near  40°-50°S  of  the  troposphere 
propagate  upward  and  slightly  southward  to  reach  near 
50°-60°S  in  the  mesosphere.  The  characteristic  propaga¬ 
tion  results  in  the  westward  forcing  having  a  peak  near 
50°-60°S  that  maintains  the  characteristic  double-peak 
structure  of  the  mesospheric  westerlies  (Sato  et  al.,  2009a; 
Watanabe  et  al.,  2009).  Similarly,  eastward  propagating 
gravity  waves  originating  near  10°-20°N  of  the  tropo¬ 
sphere  propagate  upward  and  slightly  northward  to  reach 
~  15°-25°N  of  the  mesosphere,  causing  a  eastward  forc¬ 
ing  maximum  at  these  latitudes.  The  eastward  forcing  with 
characteristic  structure  is  likely  responsible  for  the  sum¬ 
mertime  mesospheric  easterlies  tilting  toward  higher  lati¬ 
tudes  (Sato  et  al.,  2009a).  It  should  be  noted  that  gravity 


wave  forcing  is  also  dominant  in  the  summer  stratosphere. 
Watanabe  (2008)  compared  the  meridional  distributions  of 
zonal  momentum  flux  and  zonal  gravity  wave  forcing  in 
the  GWR-GCM  simulation  to  those  calculated  using  the 
Hines  gravity  wave  parametrization  scheme.  They  showed 
that  neglect  of  the  lateral  propagation  of  gravity  waves, 
including  the  effects  of  lateral  spreading  of  the  gravity 
wave  momentum  fluxes  from  the  source  locations  in  the 
gravity  wave  drag  parametrization,  could  cause  significant 
biases  in  the  parametrized  gravity  wave  forcing. 

4.3  Seasonal  variation  of  gravity  waves  and  their  sources 

An  advantage  of  global  model  studies  is  the  ability  to 
identify  gravity  wave  sources  in  the  model  (Sato  et  al., 
2009a).  The  sources  of  gravity  waves  that  act  to  maintain 
the  weak  wind  layer  in  the  upper  mesosphere  have  not 
been  known  up  until  now. 

It  is  known  from  radar  observations  in  middle  lati¬ 
tudes  that  gravity  wave  activity  in  the  mesosphere  exhibits 
a  semiannual  variation  having  two  maxima  at  solstices, 
while  that  in  the  lower  stratosphere  shows  an  annual  vari¬ 
ation  with  a  winter  maximum  (Vincent  and  Reid  1983: 
Tsuda  et  al.,  1990;  Sato  1994).  Observations  also  indi¬ 
cate  that  the  momentum  flux  pu'w'  in  the  mesosphere 
is  positive  in  summer  and  negative  in  winter.  A  likely 
explanation  for  the  seasonal  cycle  is  critical  level  filter¬ 
ing  of  gravity  waves  in  the  seasonally  varying  background 
winds  (e.g.,  Fukao,  2007).  The  seasonal  cycle  of  momen¬ 
tum  fluxes  simulated  by  the  high-resolution  model  (GWR- 
GCM,  Watanabe  et  al.,  2008)  is  shown  as  a  function  of 
latitude  in  the  mesosphere  and  lower  stratosphere  in  Fig¬ 
ure  10  (Sato  et  al.,  2009a).  It  is  seen  that  pu'w'  exhibits  an 
annual  variation  which  is  negative  in  winter  and  positive 
in  summer  both  in  the  lower  stratosphere  and  the  meso¬ 
sphere,  and  the  seasonal  cycle  in  mid-latitudes  is  con¬ 
sistent  with  radar  observations.  In  the  lower  stratosphere, 
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Figure  9.  E-P  flux  vectors  (arrows)  and  eastward  accelerations  of  zonal  mean  zonal  wind  due  to  divergence  of  E-P  flux  (colors)  for  July 
(average),  (a)  Total  wave  components,  (b)  PW  group  with  zonal  wavenumber  (fc)  =1-3,  (c)  MW  group  with  horizontal  wavelengths  greater 
than  950  km  excluding  the  k  =  1-3  components,  (d)  GW  group  with  horizontal  wavelengths  of  188-930  km.  Vertical  component  of  E-P  flux 
is  multiplied  by  250.  Scales  of  arrows  are  modified  for  clarity,  depending  on  magnitude  of  E-P  flux.  Color  scale  is  logarithmic.  Contours 
denote  zonal  mean  zonal  wind  in  10  m  s_1  intervals  (Watanabe  et  al.,  2008). 


flux  maxima  move  from  regions  with  high  topography  in 
winter  to  monsoon  convective  regions  in  summer,  sug¬ 
gesting  different  sources  in  winter  and  in  summer.  An 
interesting  feature  is  that  the  latitudes  of  the  momentum 
flux  maxima  at  the  respective  solstices  are  higher  in  the 
mesosphere  than  in  the  lower  stratosphere.  This  is  partly 
due  to  gravity  wave  filtering  by  the  mean  wind,  but  also 
due  to  gravity  wave  propagation  toward  the  jet  core,  as 
shown  in  the  next  section.  Thus,  it  is  found  from  the  model 
results  that  the  annual  variation  of  gravity  wave  fluxes  in 


the  mesosphere  is  caused  by  both  seasonal  source  varia¬ 
tions  and  by  critical  level  filtering. 

The  horizontal  distribution  of  gravity  wave  sources 
is  examined  using  a  horizontal  map  of  pu'w'.  Results 
for  July  are  shown  in  Fig.  11.  In  the  winter  (southern) 
hemisphere,  there  are  several  regions  with  largely  nega¬ 
tive  pu'w'  around  the  high  mountains  such  as  the  Andes, 
Antarctic  Peninsula,  and  the  east  coast  of  Australia.  In 
addition  to  these  isolated  regions,  negative  pu'w'  is  dis¬ 
tributed  zonally  along  the  eastward  jet.  These  may  be 
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Figure  11.  A  horizontal  map  of  momentum  flux  associated  with 
gravity  waves  (colours)  at  100  hPa  in  the  lower  stratosphere.  Red 
solid,  red  dashed  and  black  contours  indicate  zero,  westward  and 
eastward  winds  in  m  s-1  (Sato  et  al.,  2009a). 


assumptions  used  in  most  gravity  wave  parametrizations. 
If  these  model  variations  in  sources  are  realistic,  they  may 
be  important  in  explaining  the  asymmetric  structure  of  the 


mesospheric  jets  in  winter  and  summer,  as  described  in 

Figure  10.  Time-latitude  cross-sections  of  momentum  flux  associ-  section  4  4 
ated  with  gravity  waves  (colours)  and  longitudinal  mean  eastward 
winds  (contours)  at  0.03  hPa  in  the  upper  mesosphere  and  at  100 
hPa  in  the  lower  stratosphere  (Sato  et  al.,  2009a) 

4.4  Three  dimensional  propagation  of  gravity  waves 


due  to  the  gravity  waves  excited  by  spontaneous  adjust¬ 
ment  around  jet-frontal  systems  (OSullivan  and  Dunker- 
ton,  1995;  Kawatani  et  al.  2004;  Zhang,  2004;  Plougonven 
and  Snyder.  2007;  Sugimoto  et  al.  2008;  Sato  and  Yoshiki, 
2008;  Tateno  and  Sato,  2008).  In  the  summer  (Northern) 
hemisphere,  strong  positive  pu'w'  is  observed  only  in  the 
Indian  and  African  monsoon  regions.  Westward  winds  are 
dominant  in  the  lower  stratosphere  of  this  region.  Gravity 
waves  having  positive  pu'w'  generated  by  vigorous  con¬ 
vection  generally  penetrate  into  the  middle  atmosphere 
only  through  this  westward  wind  region,  an  open  window 
to  the  middle  atmosphere.  In  January,  a  similar  dominance 
of  positive  pu’w’  is  observed  in  the  Australian  monsoon 
region  (not  shown).  These  model  results  suggest  that  grav¬ 
ity  waves  sources  are  not  distributed  uniformly  in  the  hor¬ 
izontal  and  have  significant  seasonal  variations,  unlike  the 


As  illustrated  in  most  standard  text  books, 
topographically-forced  gravity  waves  in  a  steady 
wind  regime  have  zero  horizontal  group  velocity.  OGWD 
parametrizations  assume  this  explicitly,  however  zero 
horizontal  group  velocity  is  true  only  in  a  two  dimen¬ 
sional  hydrostatic  framework.  In  three  dimensional  space, 
topographically-forced  gravity  waves  propagate  leeward 
(e.g.  Shutts,  1995;  Preusse  et  al.  2002).  This  is  because 
the  energy  of  topographically-forced  gravity  waves  can 
be  freely  advected  by  the  background  wind  component 
perpendicular  to  the  wavenumber  vector,  although  the 
intrinsic  group  velocity  is  balanced  with  the  background 
wind  component  parallel  to  the  wavenumber  vector  (Sato 
et  al.  2009b).  The  model  simulation  and  a  ray-tracing 
analysis  indicate  that  with  this  mechanism,  gravity  waves 
from  the  Andes  propagate  leeward  over  a  longitudinal 
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distance  of  about  70°.  Another  important  aspect  of  lateral 
propagation  is  the  focusing  of  gravity  waves  into  the  jet. 
The  modification  of  the  wavenumber  vector  by  the  latitu¬ 
dinal  gradient  of  the  mean  wind  described  by  Dunkerton 
(1984)  and  Preusse  et  al.  (2002)  acts  on  gravity  waves 
with  negative  pu'w'  in  the  eastward  jet.  Similar  wave 
focusing  occurs  for  gravity  waves  with  positive  pu'w'  in 
westward  jets  (Eckermann,  1992;  Sato  et  al.,  2009a). 

4.5  Gravity  waves  as  a  driving  force  for  the  QBO 

Atmospheric  general  circulation  models  (AGCMs)  are 
effective  tools  with  which  to  study  the  roles  of  atmos¬ 
pheric  waves  in  driving  the  QBO  (Takahashi  1996,  1999; 
Horinouchi  and  Yoden  1998;  Hamilton  et  al.,  1999,  2001; 
Giorgetta  et  al.,  2002,  2006;  Shibata  and  Deushi  2005; 
Kawatani  et  al.,  2005,  2010a,  2010b).  Kawatani  et  al. 
(2010a)  examined  the  relative  contribution  of  equatori- 
ally  trapped  wave  modes  (EQWs)  and  internal  gravity 
waves  in  driving  the  QBO  using  outputs  over  3  model 
years  from  the  T213L256  GWR-GCM  (Watanabe  et  al., 
2008).  Here  EQWs  are  defined  with  order  n  =  —  1  to  2 
with  zonal  wavenumber  s  <  11  and  equivalent  depths  of 
2-90m,  while  gravity  waves  have  s  >  12.  (Kawatani  et 
al.,  2010a)  Figure  12a  presents  a  time-height  cross  sec¬ 
tion  of  the  EP-flux  divergence  due  to  all  resolved  wave 
components  averaged  for  10°S-10°N.  A  QBO-like  oscil¬ 
lation  having  a  period  of  approximately  15  months  is  obvi¬ 
ous.  The  period  is  shorter  than  the  observed  value,  prob¬ 
ably  partly  due  to  the  model's  underestimation  of  mean 
ascent  in  the  tropics.  The  amplitude  and  bottom  levels  of 
the  QBO-like  oscillation  are  realistic,  however.  It  is  also 
clear  that  eastward  wave  forcing  corresponds  well  to  the 
eastward  shear  of  the  QBO,  while  westward  forcing  cor¬ 
responds  to  westward  shear,  indicating  that  the  waves  in 
the  model  certainly  drive  the  QBO-like  oscillation.  It  is 


important  to  note  that  the  relationship  between  wave  forc¬ 
ing  and  vertical  zonal  wind  shear  is  also  clear  in  the  SAO 
around  the  stratopause. 

Figure  12b  show  the  time  series  of  the  EP-flux 
divergence  due  to  all  waves,  eastward  EQWs,  westward 
EQWs,  and  gravity  waves  as  well  as  the  forcing  due  to 
the  residual  circulation  at  30  hPa  averaged  for  10°S- 
10°N.  Generally,  the  forcing  due  to  the  residual  circulation 
will  be  opposite  to,  and  smaller  than,  the  total  wave 
forcing.  In  the  eastward  shear  phase,  peaks  of  eastward 
forcing  due  to  the  eastward  EQWs  and  gravity  waves  are 
almost  coincident.  The  contributions  of  eastward  EQWs 
and  gravity  waves  to  the  eastward  forcing  are  ~25- 
50%  and  ~50-75%,  respectively.  In  the  westward  shear 
phase,  westward  EQWs  contribute  ^10%  at  most  to  QBO 
driving.  On  the  other  hand,  the  contribution  by  Rossby 
waves  propagating  from  the  winter  hemisphere  was  about 
10-20%  (not  shown).  Consequently,  gravity  waves  play 
a  crucial  role  in  driving  the  QBO  in  its  westward  shear 
phase. 

In  order  to  investigate  the  horizontal  scales  of  gravity 
waves  driving  the  QBO,  time  series  of  EP-flux  divergence 
are  shown  separately  for  the  components  of  12  <  s  <  42, 
43  <  s  <  106,  and  107  <  s  <  213  (Fig.  12c).  The  east¬ 
ward  forcing  due  to  each  component  is  comparable  in  the 
eastward  shear  phase.  In  contrast,  components  with  42  < 
s  <  213  (horizontal  wavelength  <~1000  km)  are  dom¬ 
inant  for  westward  forcing  in  the  QBO  westward  shear 
phase.  These  results  are  in  general  agreement  with  earlier 
studies  using  parametrized  gravity  waves  (Giorgetta  et  al., 
2002),  but  provide  far  more  detail. 

Recent  satellite  and  modelling  studies  indicate  that 
wave  activity  depends  greatly  on  both  latitude  and  longi¬ 
tude  in  the  equatorial  region  (Alexander  et  al.,  2008;  Ern 
et  al.,  2008;  Kawatani  et  al.,  2009).  Using  outputs  from 
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Zonal  mean  zonal  wind  &  EP-flux  div  1  OS- ION 


to  inertia-gravity  waves  derived  by  Miyahara  (2006).  It  is 
shown  that  the  Walker  circulation  plays  a  crucial  role  in 
filtering  waves  propagating  into  the  lower  stratosphere  - 
westward  winds  in  the  Eastern  Hemisphere  of  the  upper 
troposphere  prevent  westward  waves  from  entering  the 
stratosphere,  and  eastward  winds  in  the  Western  Hemi- 


Figure  12.  Time-height  cross  section  of  zonal-mean  zonal  wind 
(contour)  and  EP-flux  divergence  (shaded)  at  10°S10°N.  Red  and 
blue  colors  correspond  to  eastward  and  westward  forcing,  respec¬ 
tively.  The  contour  and  shading  intervals  are  5  m  s"1  and  0.1  m 
s'1  day-1,  respectively.  Eastward  and  westward  winds  are  shown 
with  solid  and  dashed  lines,  (b)  Time  variation  of  EP-flux  diver¬ 
gence  due  to  all  waves  (black),  eastward  EQWs  (blue),  westward 
EQWs  (green),  gravity  waves  (red),  and  forcing  due  to  the  resid¬ 
ual  circulation  (yellow)  at  30  hPa  averaged  front  10°S  to  10°S.  (c) 
The  same  as  (b).  but  for  EP-flux  divergence  due  to  12  <  s  <  42 
(green),  43  <  s  <  106  (yellow),  and  107  <  s  <  213  (red).  Note 
that  the  range  of  the  ordinate  axis  of  (c)  is  different  from  that  of  (b). 

(After  Kawatani  et  al.  (2010a)) 


sphere  block  eastward  waves  -  a  generalization  of  the 
Kelvin  wave  result  reported  by  Ryu  et  al.  (2008).  In  the 
eastward  shear  zone  of  the  QBO,  the  eastward  forcing  due 
to  gravity  waves  in  the  Eastern  Hemisphere  is  much  larger 
than  that  in  the  Western  Hemisphere.  On  the  other  hand, 
in  the  westward  wind  shear  zone,  westward  wave  forcing 
does  not  vary  much  in  the  zonal  direction,  despite  slightly 
larger  westward  wave  forcing  in  the  Western  Hemisphere 
than  in  the  Eastern  Hemisphere.  Kawatani  et  al.  (2010b) 
concluded  that  the  zonal  variation  of  wave  forcing  in  the 
stratosphere  results  from  three  factors:  (1)  zonal  variation 
of  wave  sources,  (2)  vertically  sheared  zonal  winds  asso¬ 
ciated  with  the  Walker  circulation,  and  (3)  the  phase  of  the 
QBO. 

5  Data  assimilation  techniques  for  estimating  the 
Gravity  wave  momentum  budget 


the  T213L256  AGCM  Kawatani  et  al.  (2010b)  focused  on 
the  3D  distribution  of  wave  forcing.  In  both  eastward  and 
westward  shear  phases  of  the  QBO,  EP-flux  divergence 
due  to  gravity  waves  is  due  mainly  to  the  vertical  compo¬ 
nent  of  EP-flux,  which  suggests  that  a  gravity  wave  drag 
parametrization  assuming  only  vertical  wave  propagation 
is  appropriate,  at  least  for  a  QBO  simulation.  However, 
the  divergence  of  meridional  components  of  EP-flux  asso¬ 
ciated  with  EQWs  is  comparable  to  the  contribution  of 
vertical  components  where  vertical  shear  in  QBO  zonal 
winds  is  large.  The  longitudinal  dependence  of  wave  forc¬ 
ing  was  examined  using  3D  wave  activity  flux  applicable 


In  section  3,  high  resolution  observations  were  used  to 
determine  the  momentum  flux  transported  by  small-scale 
waves.  An  alternative  approach  is  to  estimate  the  miss¬ 
ing  force  in  low-resolution  analyses  using  time-averaged 
observations  and  zonal-mean  balance  conditions  (e.g. 
Shine  (1989),  Marks  (1989))  and  resolved  wave  forc¬ 
ing  (Alexander  and  Rosenlof,  1996).  More  recently,  data 
assimilation  methods  have  been  used  to  capture  the  3D, 
time  evolving  drag  field.  The  actual  missing  force  due  to 
unresolved  waves  is  the  result  of  interactions  between  the 
mean  flow,  the  resolved  waves  and  the  unresolved  waves. 
In  particular,  the  forcing  by  unresolved  waves  changes  the 
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Figure  13.  Zonal  momentum  flux  (contour  interval  is  0.01  Pa,  positive  zonal  flux  is  shaded)  at  100  hPa  for  January  and  September  2002. 

After  Pulido  and  Thuburn  (2008). 


mean  flow  which  then  produces  changes  in  the  forcing  by 
the  resolved  waves.  By  constraining  the  large  scale  flow 
with  measurements  (using  data  assimilation),  the  forcing 
due  to  resolved  waves  is  known  so  that  the  misfit  between 
the  observed  and  forecasted  drag  can  be  attributed  to  unre¬ 
solved  waves.  Thus,  data  assimilation  methods  have  the 
unique  advantage  of  being  able  to  separate  the  response  of 
the  drag  due  to  resolved  waves  from  that  due  to  unresolved 
waves.  In  this  section,  the  standard  assimilation  problem  is 
introduced  before  considering  its  application  to  the  inverse 
problem  of  estimating  the  missing  force  due  to  unresolved 
gravity  waves. 


5. 1  The  data  assimilation  problem 

Data  assimilation  combines  measurements  and  model 
forecasts  according  to  their  accuracies  and  produces  a 
’’best”  estimate  of  the  state  of  the  atmosphere  on  a  model 
grid.  This  definition  is  appropriate  to  the  process  per¬ 
formed  at  numerical  weather  prediction  centres.  Mathe¬ 
matically,  the  process  can  be  represented  as  follows.  The 
assimilation  step  seeks  a  model  state,  x,  which  best  fits 
observations,  z,  and  which  is  not  too  far  from  a  back¬ 
ground  model  state,  Xb.  This  can  be  done  by  minimizing 


the  following  cost  function: 

J (x)  =  ~ (x  -  xb)TB_1(x  -  x6)  +  ~ (z  -  iT(x))5  R_1(z  -  fT(x)). 

(6) 

Here  H  is  a  possibly  nonlinear  operator  mapping  model 
to  observed  variables,  and  B  and  R  are  error  covariance 
matrices  corresponding  to  the  model  and  observed  state 
vectors.  The  background  model  state,  xb,  typically  comes 
from  a  short  term  (6  or  12-h)  forecast  valid  at  the  time  of 
the  analysis.  If  the  accuracy  of  measurements  and  back¬ 
ground  state  is  known,  then  the  state  which  minimizes  (6) 
is  optimal  in  that  it  minimizes  the  analysis  error  variance. 

However,  error  covariances  (particularly  those  associated 
with  the  background  state)  can  never  really  be  known 
and  must,  at  best,  be  modelled  approximately.  Thus  the 
optimality  principle  holds  only  in  theory.  Once  an  ini¬ 
tial  state  is  obtained,  a  short  term  model  forecast  can  be 
launched  to  obtain  a  background  state  for  a  later  assimi¬ 
lation.  Then  the  assimilation  cycles  continuously  every  6 
or  12  hours  through  alternating  assimilation  and  forecast 
steps.  (6)  rather  precisely  defines  the  so-called  3D  vari¬ 
ational  assimilation  procedure  (or  3D-Var).  By  allowing 
H  to  include  temporal  interpolation  of  the  model  state  to 
the  time  of  observations  distributed  over  a  time  window 
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using  the  forecast  model,  this  equation  can  also  describe 
the  so-called  4D  variational  assimilation  procedure  (or 
4D-Var).  Finally,  note  that  Optimal  Interpolation  and  3D- 
Var  are  formally  equivalent  in  the  case  of  linear  observa¬ 
tion  operators  and  Gaussian  observation  and  background 
errors.  Thus  (6)  can  also  describe  the  Optimal  Interpola¬ 
tion  scheme. 

5 .2  Evidence  that  the  mesosphere  is  constrained  by  mea¬ 
surements  from  the  troposphere  and  stratosphere 

Gravity  wave  drag  reaches  large  amplitude  in  the  upper 
mesosphere  and  lower  thermosphere  region.  Since  obser¬ 
vations  that  are  assimilated  into  weather  and  climate  mod¬ 
els  are  primarily  from  the  troposphere  and  stratosphere,  it 
is  important  to  know  whether  such  observations  can  con¬ 
strain  the  mesosphere.  Conversely,  if  the  mesosphere  is 
slaved  to  the  lower  atmosphere  then  the  misfit  of  mea¬ 
surements  and  forecasts  in  the  mesosphere  can  be  used 
to  retrieve  information  about  the  troposphere  and  strato¬ 
sphere.  Recently,  Ren  et  al.  (2008)  have  shown  that  fore¬ 
casts  of  the  2002  southern  hemisphere  stratospheric  sud¬ 
den  warming  that  capture  the  vortex  splitting  event  also 
produce  a  mesospheric  cooling  while  forecasts  that  do 
not  produce  a  vortex  splitting  have  much  weaker  meso¬ 
spheric  cooling.  In  these  simulations,  observations  were 
used  only  below  1  hPa  so  that  model  forecasts  in  the 
mesosphere  were  due  to  the  model  response  to  data  inser¬ 
tion.  Of  the  mechanisms  available  to  produce  a  meso¬ 
spheric  cooling  during  forecasts,  parametrized  GWD  was 
shown  to  be  the  dominant  one.  Thus,  in  a  dramatic  event 
such  as  a  sudden  warming,  GWD  is  playing  an  important 
role  in  determining  the  mesospheric  forecast.  While  this 
results  demonstrates  that  the  mesosphere  is  slaved  to  the 
lower  atmosphere  in  at  least  one  warming,  a  more  statis¬ 
tical  result  was  obtained  by  Nezlin  et  al.  (2009).  In  this 
work,  observations  from  the  troposphere  and  stratosphere 


were  able  to  define  the  large  scales  (below  wavenumber 
10)  in  the  mesosphere.  Since  observations  were  simulated 
from  a  model  run,  there  was  an  underlying  assumption 
that  the  model  is  realistic.  Even  with  this  assumption,  and 
with  perfect  observations  (no  random  errors  added),  meso¬ 
spheric  scales  higher  than  wavenumber  10  could  not  be 
determined.  Thus  only  the  large  scale  mesospheric  flow 
can  be  constrained  by  observations  from  below  the  meso¬ 
sphere. 

The  role  of  initial  conditions  in  the  troposphere  and 
stratosphere  provided  by  data  assimilation  in  control¬ 
ling  a  subsequent  mesospheric  forecast  without  a  local 
data  assimilation  constraint  was  similarly  demonstrated  in 
studies  by  Coy  et  al.  (2005)  and  Siskind  et  al.  (2007). 
These  earlier  studies  also  validated  their  mesospheric 
forecasts  against  independent  SABER  temperature  obser¬ 
vations.  The  observational  constraints  showed  that  the 
mesospheric  response  to  stratospheric  warmings  can  be 
quite  variable  from  event  to  event,  but  the  response  is 
typically  characterized  by  a  shallow  lower  mesospheric 
cooling  layer  that  then  transitions  to  an  upper  mesospheric 
or  lower  thermospheric  secondary  warm  layer  response. 
These  details  differ  from  the  Ren  et  al.  (2008)  study,  and 
they  may  be  sensitive  to  the  details  of  the  GWD  parame¬ 
terization  used  in  the  different  models. 

5.3  Using  data  assimilation  to  estimate  gravity  wave 
drag  morphology 

The  response  of  the  large-scale  extratropical  flow  to  an 
external  localized  forcing  is  non-local  and  may  include 
geostrophic  and  transient  components  (see  Fig.  6  in 
Flaynes  et  al.  1991,  also  Fig.  1  in  Pulido  and  Thuburn 
2005).  If  the  initial  state  is  known,  then  an  initial  uncer¬ 
tain  forcing  added  to  the  momentum  equations  will  lead 
to  a  departure  of  the  predicted  drag  (or  wind)  from  the 
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observed  drag  (or  wind).  Inverse  techniques  use  this  mis¬ 
match  of  the  observed  and  predicted  drag  and  knowledge 
of  the  flow  evolution  (model  equations)  to  trace  back  and 
recover  the  initial  localized  forcing.  Pulido  and  Thuburn 
(2005)  were  the  first  to  apply  an  inverse  technique  based 
on  the  concepts  of  4D  variational  data  assimilation  to  esti¬ 
mate  GWD.  A  cost  function  that  measures  the  difference 
between  the  state  of  the  dynamical  model  and  observa¬ 
tions  was  defined: 

J(x0,X)  =  ^(zi  -  Mi (x0 , X) )  TR+  1  (z-i  -  M,(x0,X)) 

i 

(7) 

where  z,  are  the  observations  (the  analysis  in  this  context) 
at  time  i,  x0  is  the  initial  condition,  X  is  the  GWD,  and  M 
is  the  dynamical  model  so  that  M,  (x0 ,  X)  gives  the  state 
of  the  model  at  time  i.  H  does  not  appear  explicitly  in  (7) 
because  observations  z,  were  first  transformed  to  model 
state  variables.  The  assumption  of  complete  ignorance 
of  the  gravity  wave  drag  is  taken  so  that  there  is  no 
background  term  (first  term  in  (7))  in  the  cost  function.  As 
the  observations  and  the  initial  conditions  x0  are  known, 
the  only  unknown  field  in  the  cost  function  (7)  is  the  GWD 
X.  The  determination  of  the  minimum  of  the  cost  function 
gives  the  optimum  gravity  wave  drag. 

Budget  techniques  give  a  gravity  wave  deceleration 
centre  in  the  winter  hemisphere  and  a  weaker  deceleration 
centre  in  the  summer  hemisphere  above  the  stratopause, 
while  gravity  waves  accelerate  the  jets  in  the  stratosphere 
(e.g.  Alexander  and  Rosenlof,  1996).  With  data  assim¬ 
ilation,  the  estimated  zonal  mean  of  the  zonal  compo¬ 
nent  of  GWD  qualitatively  resembles  the  results  obtained 
by  budget  studies.  However,  there  are  quantitative  differ¬ 
ences  such  as  stronger  and  more  concentrated  deceleration 
centres  being  found  with  the  data  assimilation  technique 


(Pulido  and  Thuburn  2006).  While  the  source  of  these  dif¬ 
ferences  remains  unclear,  one  advantage  of  assimilation 
techniques  over  budget  studies  is  clear:  that  they  are  able 
to  infer  3D,  day-to-day  GWD  fields,  of  both  zonal  and 
meridional  components. 

Using  the  GWD  estimated  with  data  assimilation, 
Pulido  and  Thuburn  (2008)  inferred  the  momentum  flux 
at  100  hPa  assuming  that  momentum  fluxes  are  negli¬ 
gible  above  0.24  hPa.  The  calculated  momentum  flux  is 
dominated  by  the  GWD  in  the  lower  stratosphere  where 
model  errors  may  have  some  impact  on  the  GWD  esti¬ 
mation.  Figure  13  shows  the  monthly  mean  momentum 
flux  at  100  hPa  for  September  2002.  Maximum  momen¬ 
tum  fluxes  are  found  over  Antarctica  (~  +50  mPa)  and 
over  South  America  (~  —20  mPa)  which  could  be  related 
to  the  orography  or  the  storm  track  location.  During  boreal 
winter,  maximum  momentum  fluxes  are  located  over  Asia 
(~  —40  mPa). 

The  assimilation-estimated  GWD  also  has  reasonable 
features  at  low  latitudes  where  GWD  is  expected  to 
play  a  role  in  the  quasi-biennial  oscillation  (QBO).  With 
a  one-year  GWD  estimate,  Pulido  and  Thuburn  (2008) 
noted  that  at  low  latitudes  the  forcing  changes  from 
positive  to  negative  at  40  hPa  (both  are  accelerating  the 
magnitude  of  zonal  mean  flow,  see  Fig.  14)  in  consonance 
with  the  zonal  wind  that  is  changing  from  the  positive 
to  the  negative  phase  of  the  QBO.  At  higher  altitudes 
a  semiannual  oscillation  is  found  in  the  zonal  GWD, 
however  the  GWD  has  opposite  phase  to  the  zonal  wind 
semiannual  oscillation,  and  is  superimposed  on  a  mean 
positive  forcing  (left  and  middle  panels  of  Fig.  14). 

Estimation  of  GWD  with  data  assimilation  relies  on 
two  major  assumptions:  The  forcing  produced  by  gravity 
waves  must  be  larger  than  model  error  and  it  must  be  a 
large-scale  and  systematic  forcing  so  that  low-resolution 
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Figure  14.  Daily  zonal  mean  GWD  (continuous  line)  estimated  at  the  Equator  with  a  data  assimilation  technique  for  2002,  zonal  mean 
wind  (dotted  line)  and  linear  regression  of  the  drag  (dashed  line)  at  0.24  hPa  (left),  1.15  hPa  (middle)  and  40  hPa  (right  panel).  After  Pulido 

and  Thuburn,  2008. 


observations  can  capture  the  response  of  the  flow  to  the 
forcing.  The  latter  assumption  has  a  physical  foundation: 
the  systematic  filtering  of  gravity  waves  by  background 
winds  in  the  lower  atmosphere  produces  a  large  scale, 
slow  forcing  in  the  mesosphere.  The  former  assumption 
is  justifiable  since  climate  models  without  a  gravity  wave 
parametrization  or  Rayleigh  friction  are  expected  to  con¬ 
tain  a  strong  bias  above  the  stratopause,  the  so-called 
missing  force  (e.g.  Holton  1982).  Since  a  poor  representa¬ 
tion  of  other  physical  processes  in  models  (e.g.  radiation) 
will  also  be  attributed  to  GWD,  it  is  important  that  GWD 
be  the  dominant  source  of  model  error.  The  phasing  of 
these  annual  and  semiannual  variations  are  in  agreement 
with  monthly  means  derived  in  Alexander  and  Rosenlof 
(2003). 

Thus  far,  assimilation  techniques  have  been  applied 
to  estimate  the  drag  due  to  unresolved  gravity  waves.  In 
the  future,  these  techniques  may  also  be  applied  to  directly 
determine  parameters  used  in  GWD  schemes.  The  use  of 
data  assimilation  or  inverse  techniques  offers  an  objective 
way  to  estimate  climate  model  parameters.  Such  an  objec¬ 
tive  parameter  estimation  may  be  of  paramount  impor¬ 
tance  to  realistically  represent  the  QBO,  a  circulation  fea¬ 
ture  that  is  still  not  well  represented  in  state-of-the-art 
climate  models. 
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6  Summary 

Here  we  summarize  a  few  key  points  regarding  the  effects 
of  gravity  waves  in  climate  models,  results  from  obser¬ 
vational  analyses  that  constrain  these  effects,  and  very 
high-resolution  model  studies  that  highlight  both  the  capa¬ 
bilities  of  these  models  as  well  as  the  challenges  that  still 
lie  ahead  for  realistic  treatment  of  gravity  wave  effects  on 
the  large-scale  circulation. 

Parametrized  orographic  gravity  wave  drag  (OGWD) 
has  strong  effects  on  mean  winds  and  planetary  wave 
propagation  in  the  winter  hemisphere  stratosphere,  and 
it  is  through  this  mechanism  that  OGWD  affects  climate 
change  response  patterns  in  models.  OGWD  also  has  a 
strong  effect  on  variability  both  in  the  stratosphere  and 
at  the  surface.  An  important  tuning  parameter  controlling 
these  effects  (e),  which  scales  the  drag  profile,  remains  dif¬ 
ficult  to  quantify  with  observations.  Angular  momentum 
conservation  is  a  critical  constraint  that  must  be  applied 
to  ensure  robust  results  in  climate  change  studies  that  use 
parametrized  GWD. 

Recent  studies  have  noted  parametrized  GWD  effects 
on  some  climate  responses  in  models.  OGWD  affects 
predicted  trends  in  the  Brewer-Dobson  circulation  in  cli¬ 
mate  change  scenarios.  Changes  in  this  circulation  have 
repercussions  for  the  transport  of  trace  constituents  and 
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for  ozone  chemistry.  Parametrized  non-orographic  gravity 
wave  drag  (NGWD)  plays  a  crucial  role  in  chemistry- 
climate  forecasts  of  ozone  recovery  where  it  is  used  to 
control  polar  temperatures  in  model  studies,  particularly 
in  the  southern  hemisphere.  Small  changes  to  the  tuning 
parameter  that  sets  the  source-level  momentum  flux  can 
have  very  large  effects  on  the  seasonal  development  of  the 
ozone  hole.  Source  parametrizations  for  known  sources 
such  as  convection,  fronts,  and  jet  imbalance  are  needed 
to  give  gravity  waves  sensitivity  to  meteorology  in  cli¬ 
mate  models  so  they  can  evolve  with  changing  climate. 
Although  experimental  versions  of  such  parametrizations 
exist,  they  are  currently  poorly  constrained. 

Global-scale  disUibutions  of  gravity  wave  momen¬ 
tum  flux  and  other  propagation  properties  have  been 
derived  from  both  satellite-  and  balloon-borne  measure¬ 
ments.  Each  measurement  and  analysis  method  observes 
only  a  portion  of  the  full  spectrum  of  momentum  fluxes, 
but  the  combined  set  of  measurements  can  provide  very 
good  spectral  coverage.  Momentum  fluxes  derived  from 
observations  can  also  depend  on  the  resolution  at  which 
they  are  reported  because  wave  events  can  be  very  local¬ 
ized  in  both  space  and  time,  a  concept  termed  intermit- 
tency.  Because  of  intermittency,  local  values  can  be  more 
than  an  order  of  magnitude  larger  than  areal  or  time- 
averaged  fluxes. 

Time-averaged  global  observations  of  momentum 
fluxes  derived  from  different  measurements  are  showing 
some  similarities:  Largest  fluxes  in  the  winter  hemisphere 
over  topography,  enhanced  fluxes  at  low  latitudes  and 
in  the  summer  hemisphere  over  known  centers  of  deep 
convection.  Uncertainties  in  magnitude  and  direction  of 
propagation  remain,  however  methods  for  intercompari¬ 
son  and  correction  for  various  limitations  of  the  differ¬ 
ent  observations  are  being  developed  that  should  allow 


a  more  quantitative  global  picture  to  emerge  in  the  very 
near  future.  Future  measurements  with  improved  resolu¬ 
tion  that  extend  the  global  coverage  of  the  balloon  data 
and  that  quantify  the  wave  propagation  directions  in  satel¬ 
lite  data  will  have  major  impacts  on  reducing  uncertain¬ 
ties. 

Satellite  evidence  suggests  OGWD  may  play  a  larger 
role  in  the  southern  hemisphere  than  previously  thought 
because  of  small  island  sources  of  mountain  waves.  The 
Southern  Andes  and  Antarctic  Peninsula  also  appear  in 
all  the  observations  and  models  as  a  locus  of  maximum 
fluxes.  Nonhydrostatic  propagation  effects  and  nonlinear 
processes,  both  of  which  are  neglected  in  current  OGWD 
parametrizations,  likely  play  important  roles  here.  Future 
work  may  lead  to  better  methods  for  treatment  of  these 
processes  in  the  parametrizations. 

The  most  recent  gravity-wave  resolving  high- 
resolution  global  climate  models  have  achieved  a  high 
level  of  realism,  and  the  gravity  wave  momentum  fluxes 
in  these  models  show  excellent  agreement  with  obser¬ 
vations  from  radar  and  satellite.  This  is  somewhat  sur¬ 
prising  because  observational  evidence  suggests  that  a 
large  fraction  of  the  gravity  wave  momentum  flux  remains 
unresolved  in  these  models.  The  3D  picture  afforded  by 
the  model  allows  calculation  of  global  distributions  of 
the  vector  momentum  flux,  the  associated  GWD,  and  the 
sources  of  the  gravity  waves.  The  high-resolution  model 
results  suggest  that  waves  in  the  gravity  wave  part  of  the 
spectrum  provide  more  than  half  of  the  wave  forcing  nec¬ 
essary  to  drive  the  QBO. 

Patterns  in  the  gravity  wave  momentum  flux  derived 
from  data  assimilation  methods  show  some  consistency 
with  the  observations  and  similarity  to  patterns  derived 
from  high  resolution  models.  These  similarities  are 
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encouraging  since  the  methods  are  all  dramatically  dif¬ 
ferent.  The  similarities  among  the  global  distributions 
derived  from  these  very  different  methods  suggest  that  a 
quantitative  picture  of  the  sources  of  gravity  waves  and 
the  global  disUibutions  of  momentum  flux  will  emerge  in 
the  near  future. 

Tightening  constraints  on  a  GWD  parametrization 
will  require  the  additional  step  of  translating  these  find¬ 
ings  into  practical  limits  on  the  set  of  tuning  parameters, 
which  differs  from  one  parameterization  to  another.  This 
can  only  be  accomplished  through  collaborative  research 
between  model  developers  and  observational  scientists. 
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